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Abstract We discuss how recent advances in observations, theory and numerical sim¬ 
ulations have allowed the stellar community to progress in its understanding of stellar 
convection, rotation and magnetism and to assess the degree to which the Sun and 
other stars share similar dynamical properties. Ensemble asteroseismology has become 
a reality with the advent of large time domain studies, especially from space missions. 
This new capability has provided improved constraints on stellar rotation and activity, 
over and above that obtained via traditional techniques such as spectropolarimetry or 
Call H&K observations. New data and surveys covering large mass and age ranges have 
provided a wide parameter space to confront theories of stellar magnetism. These new 
empirical databases are complemented by theoretical advances and improved multi- 
D simulations of stellar dynamos. We trace these pathways through which a lucid 
and more detailed picture of magnetohydrodynamics of solar-like stars is beginning to 
emerge and discuss future prospects. 
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1 Solar Dynamics and Magnetism in a Stellar Context 


Time-dependent stellar magnetic fields, such as that of the Sun, are generated in their 
interior through a magnetohydrodynamic (MHD) dyna mo mechanism that involves 
inter actions between plasma flows and magnetic fields le.g. iMoffatdllQTil i lCharbonneaul 
boifll 'I. The generated magnetic fields are manifested at the surface through starspots or 
polar caps, and contribute to the dynamic activity of the star. This activity governs the 
structure of stellar co ronae, the stellar radiative flux ou tput, and transient but energetic 
events such as flares I GudelllioOdI : ISchriiver et al.ll2012l l. Stellar activity can also have a 


strong impact on the properties of stellar astrospheres which may host planets ( WoodI 


I 2 OO 4 I I. Magnetized solar-like winds extract angular momentum which drives angular 
momentum in low mass stars. Understanding the origin of stellar magnetism is therefore 
of crucial importance, and remains an outstanding problem in astrophysics. 

The Sun, our nearest star, provides easy access to spatially resolved observations 
of magnetic activity and thus sets constraints on the dynamo origin of stellar mag¬ 
netic fields. However, significant uncertainties in our current understanding of the so¬ 
lar dynamo mechanism motivates us to confront our theoretical ideas of magnetic 


vations dNoves et al. 

1984a: Ba 

unas et al.l 1995: Wright et a 

.| 2004 

1: GiampaDall2005l: 

iBohm-Vitensd 20071: 

Olah et al. 

2 OO 9 I: Donati and Landstreetl 


R,einersll2012l. and 


references therein). While, in principle, different stars with different properties may 
host distinct classes of MHD dynamos, a good starting assumption is that a reasonable 
model of solar and stellar magnetism should be able to reproduce the range of stellar 
activity taking into account appropriate changes in the internal driving mechanisms. 
For example, high mass stars have shallow surface convective zones with high turbu¬ 
lent fluxes relative to low mass stars. Stars with different rotation rates and ages will 
have different magneto-convective properties (usually characterized by varying non- 
dimensional parameters such as the Reynolds or Rossby numbers). To reduce this 
problem to a practical exercise, it is also useful to deal with solar-like stars (narrowing 
the range of spectral type consi dered), within which, there is substant ial variation in 
the nature of magnetic output ( Petit et al.ll20oA iMarsden et allbol.lh . The observa¬ 


tion and characterization of stellar convection, rotation and magnetism has undergone 
a recent boom with the outcome of better ground-based instruments and space-born 
missions. For instance, until recently chromospheric proxies such as emission in the 
core of Call H & K lines and white light modu lation were the most commonly use d 


techniques for characterizing stellar magnetism ( Wright et al. I |2004 lGiampa^l2005l i . 


The development of very sensitive spectropolarimeter such as ESPaDOnS or NARVAL 
has helped to further improve our understanding of stellar magnetism by providing 
detailed magnetic maps and sometimes the surface diffe rential rotation of the star 
( Donati et al.ll 199'^ : Ijardine et al.lll999l : IPonati et al.ll200^ . With the success of planet 


transit/asteroseis mology missions su ch as Corot, Kepler, the picture can now be re¬ 
fined even further ( Garcia et al.l[20l3 i. We have now entered an era in which systematic 


coupled analyses of chromospheric activity proxies, spectropolarimetry maps and as- 
teroseismic data is being performed. This allows us to better characterize magnetism of 


solar-like stars and its variability on short and relat i vely long terms (IBe rdvuginal 


2005 


Peti^et_al ] l2008l : iMorgenthaler et al.l l2011alf2011bl : iMorin et alJ l201,^ : iMathur et al.l 


201 flai l. The purpose of this review paper is to recap recent progress made, with the 


aid of asteroseismology, theory and numerical simulations, in our understanding of 
stellar dynamics and magnetism in solar-like stars. Other proxies, such as the time 
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evolution of surface lithium abundances, can also be useful to constrain stellar dy¬ 
nami c s but they are outside the s c ope of this review (for discussions see Brm et aP 


199^ Talon and Charbonnell l2004l : ISbordone et alJ I 2 OI 0 I : IDo Nascimento et alJ 201tt 


Spite et al. 2012I . and references therein). We begin by discussing recent observational 


findings with a special emphasis on ensemble asteroseismology and on stellar dynamics, 
magnetism and cyclic activity. We then present the current theories of stellar rotational 
history based on 1-D models and of stellar magnetism with an introduction to dynamo 
theory based on 2-D kinematic dynamo models and 3-D nonlinear turbulent numerical 
simulations of stellar dynamo. We continue with a discussion on how theoretical mod- 
els hold up against observations such as s caling laws and rotation-activity relationship 
( Pizzolato et aI.l[2bo3l : IWright et al.ll201ll ~) and provide a critical assessment of the cur¬ 
rent understanding. Finally we conclude by proposing possible future developments in 
asteroseismology and stellar dynamics to better constrain solar-like star rotation and 
activity. 


1.1 New Asteroseismic Diagnostics of Global Stellar Convection, Rotation and 
Magnetism 


Oscillations are intrinsically damped and stocha stically excited by the turbu l ent motion 


in th e outer conyect i ye layers of cool stars (e.g. [Goldreid^^n^Keele^lll^^j_ Balmforpil 
I 1992 I : iHoudek et al.l 11999 : IChaplin et al.l 1200^ : Belkacem et al. 20081 : Houdek 2 OI 0 I I. 


We refer to stars exhibiting such an oscillation pattern as solar-type stars hereafter. 
These oscillations are a l so present with higher amplitudes in red-giant stars (e.g. 

iHoudek and Goughl2002l : Ide Bidder et al.l[2009l : iBedding et al.l201^ : Mosser_et^l j201(]| l 


(and eyen in mores massiye B-type stars as recently shown by Belkacem et al.l 2009l i. 


The frequency of maximum power scales as {GM/ fe.g.. lKieldsen and Beddingl 
I 1995 II : giyen the wide dynamic range in stellar surface grayities the releyant periods 
span the range of minutes to weeks. 

The magnetic activity cycles in the Sun and other solar-like stars are often the 
consequence of dynamo processes induced by the interaction of three fundamental 
mechanisms: convection, rotation, and magnetic fields (see detailed discussion in §2). 
Asteroseismic observations place powerful observational constraints on each of these 
three processes as we will see in §1.4, 1.5, and 1.6. For example, activity-cycle induced 
variations of pulsation properties (mode linew idth and ampl i tudes ) in t he Sun and in 
solar- like stars were discussed, for example, bv iHoudek et al.l 1 200lll and IChaplin et ahl 
( 2007lh 

The long and uninterrupted high-precision photometric observations required for 
asteroseismology are ideal for characterizing long-term variations such as the ones 
induced by starspots crossing the visible disk of stars. Indeed, the dimming pro¬ 
duced by these spots leads to a measurable modulation in the photom etry with a 
period corresponding to the surface rotation at the active latitudes (e.g. BasrLet_alJ 
20n Garcfa^et_^ 20R Garci^et^ 20M;_ Lanz^ 20R: Walkowicz and Bas^ 201^ 


NiSareT'^ 2o5 ReiMroid^T'^ 20131 : McQuillan et al I 2 OI 4 I I in which the spots de¬ 

velop. One must keep in mind that on the Sun faculae and the small-scale field around 
over-compen s ate this dimm ing effect, leading to a net brightening at solar maximum 
( Spruitir2000l : [hMihlichll2012l f . Moreover their relative amplitude depends on how close 
to the limb the starspot is and to the stellar inclination angle with respect to the line of 
sight. Mapping the relative ratio of dimming through star spot vs. brightening through 
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faculae and small-scale field changes as a function of stellar spectral type, age, rotation 
rate and metallicity is an emerging field, and as of this paper the overall pattern is not 
firmly establish ed. Some studies ha ve started to address the physical origins of stellar 
variability (e.g. IShaoiro et ^l2014l i and found that a transition from net brightening 
to net dimming should occur around an index of ch romospheric activity logffl^Tf) of 
-4.9, a value which is consistent with observations ( Lockwood et al.l l2007l : iHall et al] 
l2009li . 

Ground based stellar variability studies have a long history, and they have tradi¬ 
tionally focused on young and active stars. The Convection Rotation and planetary 
Transits (CoRoT) and Kepler space missions have contributed continuous and precise 
data series for large samples, which has had a dramatic impact on the quantity and 
quality of stellar rotation data. CoROT and Kepler h ave extended our understanding to 
large numbers of older and less active field stars fe.g. iMosser et aDl2nn9l : [Mathur et al.l 
l2010bl . l2011bl : ICampante et al.ll201lli . Asteroseismology has also extended our under¬ 
standing of longer term cycles. In the solar case the number of spots in the Sun changes 
during the solar magnetic cycle l Hathaw^l20in(i . introducing a modulation of the 
variance in the photometric o bservations that can be used as a proxy of the mag- 
netic activity cycle in the Sun j Mathur et al.l[2009ll and in other solar-like stars (e.g. 
iGarcfa et al. I 2 OI 0 I : IChanlin et al. 2011a). As the temperature rises from the photo¬ 
sphere to the chromosphere, the c ores of very strong lines, such as Ca II H (396.8 nm ) 
and K (393.4 nm) emission lines ( Leightonlll959l : [wilsoniriQTSal : iBaliunas et al.lll9^ 1. 
develop emission features whose strength depends on the magnitude of the temperature 
inversion. However, seismic techniques can go further by probing structural and dy¬ 
namical changes beneath the photosphere as the magnetic cycles develop. Indeed, since 
the dawn of helioseismology, several auth ors noticed that the frequency of the acoustic 
modes increase n e ar the solar maximum (|WoodMd_and_^iqYgjll9^ Fossat et al.lll98^ : 
IPalle et al.l Il989l : iLibbrecht and Woodard 199fll : Anguera Gubau et al.l 1992ll . Nowa- 
days, it i s possible to measure variations in all the p-mo de properties not only in the 
Sun fe.g. IChanlin et 51l2000l : Ijimenez-Reves et al.l[2007h . but also in other stars such 
as the F star HD49933 measured by GoRoT, in which the amplitude as well as the 
frequency shifts of low-degree modes change with time ( Garcia et al.ll20ld i. 

Asteroseismic measurements can also be used to extract information about convec¬ 
tion. The p modes are observed as a power excess at distinct frequencies on top of a 
continuous Fourier spectrum. This background spectrum is likely due to t he presence of 


convective motions at all scales with granulation playing a key role (e.g. iHarvej 


1985; 


Lefebvre et aklBoOSh . but may also be due to the presence of faculae fe.g. iKaroffeTliD 


2013l i. Until now, neither classical observations nor one-dimensional analytical con¬ 


vection models could put constraints on granulation. The situation has changed with 
aster oseismic obser vations of thousand s of stars thanks to the CoRoT ( BaMinet_aD 
l2006l i and Kepler ( Borucki et ahllioidi missions. For instance. iMathur et al.l ll2011^ 
have determined the timescale and power of granulation for around 1000 Kepler red 
giants at different stages of their evolution going from the red giant branch up to the 
clump. Moreover, the depth of the convective zone -which is often the location 
used to compute a caracteristic convection overturning time e.g. Tconv ~ d^^^/vconv, 
with Uconv being a characteristic convective velocity - is a region characterized by a 
sharp change in stratification, since it changes from (super)-ad iabatic to subadiabatic, 
modulo the extent of convective penetration (e.g., IZahrJ[l99lh . Therefore, the sound 
speed undergoes an abrupt variation producing an aco ustic glitch that leaves a fin- 
gerprint in the frequencies of the acoustic modes (e.g., iGough and ThompsorJ ri988l : 
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Vorontso^J1988UGough 
Houdek and Goum 2007 


1990|) see also | Montemi-gL^ 199j^ 199^: 


IChristensen-DaisgaarT'et'lini20nr Mazumdar et al Jl2014|). 


1.2 The Link Between Helio- and Asteroseismology: Global Oscillations and 1-D 
Models 


Helioseismology has been providing a wealth of high-quality information about the 
structure and physics of the solar interior sin ce the first detection of the global five- 
minute oscillation modes i Claverie et al.lE979l L The large number of observable global 


solar oscillation modes, with spherical degrees up to several thousands, has allowed 
helioseismologists to infer with high accuracy from the measured oscillation frequency 
the solar sound speed and density as a function of solar radius. Such helioseismic inver¬ 
sion techniques make use of the fact that the displacement eige nfunctions are related 
to the adiabatic eigenfrequencie s through a variational principle i Chandrasekharj|l963l : 


iLvnden-Bell and Ostrikerill967l 'l leading to a linear relation between a small perturba¬ 
tion in the density (and squar ed sound speed) of a solar model and the corresponding 
change in the frequencies fe.g. lGoughlll97^ . Il984lj) . Before the results of the Sudbury 
Neutrino Observatory (SNO) experiment (see further down), the solar neutrino prob¬ 
lem, together with the high-quality helioseismic data collected from space, indicated 
clear problems in our understanding of either basic physics or the structure of the solar 
interior. In res ponse, there was a major susta ined effort to improve solar models; for a 
review see e.g. IChristensen-Dalsgaard ( 20041 ') and references therein. The rather good 


agreement (< 1%) between the Sun and 1-D solar models computed in the late 1990s 
( Christensen-Dalsgaard et al ] ll99fil:lBrun et al.lll998l. Il999ll was , howe ver, degraded if 


the solar photospheric abundances reported bv lAsplund et aP 1 2004l ~l were adopted. 


The p roblem still holds eve n after the re vision to a slightly higher value of metalic - 
ity by A^pluncLeL-alJ (1200911 . We refer to iGhristensen-Dalsgaard and HoudekI (l201fll l. 


IVillante'^t^n ~ l20iy and chapter 2 (e.g., iBasu et^iTl201^ 1 of the ISSI book ?? for a 
recent discussion on helioseismology constraints on metallicity/abundances. 

Still, the high-quality helioseismic data has led to many fundamental achievements 
in solar physics such as inverting the internal solar angular-velocity profile or helping to 
disentangle the long-standing solar neutrino problem. Seismic inversions have revealed 
that the non-uniform, surface differential rotation extends all the way to the base 


the radiative interior (e.g. Brown and Morrowl 

1987: Thompson 1990l: Thompson et al.l 

I 2 OO 3 I: iGarcfa et al.l I 2 OO 4 I; iMathur et al.l 20081: 

Eff-Darwich and Korzennik 2 OI 2 I 1 . The 


inferred radial interior density and sound-speed profiles also place strong constraints 
on the standard solar model and consequently on the theory of stellar structure and 
evolution. Since the neutrino fluxes predicted by solar models could not be explained 
by standard zero mass neutrino physics, particle physicists had to consider massive 
neutrinos and their associated flavor oscil lations (e.g. Tiirck-Ghieze et liDboOlIl . Thanks 
to the SNO experiment measurements ( Ahmad et al Jl2Q0l|) the scientihc community 
has since confirmed, by measuring both charged and neutral currents, that the original 
deficit of electronic neutrinos is due to their transition to another neutrino flavor as 
they travel from the solar core to Earth’s detectors. Another important issue, which 
has again become topical very recently by analyzing the high-quality data from the 
Heliospheric and Magnetic I mager (HMI) on N ASA’s Solar Dynamics Observatory 
(SDO), is the solar oblateness ( Kuhn et al.ll201^ and the consequent distortion of the 
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Sun’s gravitational field which may account, in part, for the remaining 0.2% d iscrepancy 
in the precession of Mercury’s orbit (e.g. ISchou et alJ Il998l : IGou riill2ni2ll . Detailed 
discussions about what we have learnt from helioseismology and how helioseismology 
can help t o address the remaining prob lems in sola r phys ics can be found in recent 
reviews bv IChristensen-DalsgaardI 1 20041 ') an d|G oug i] (l2013ll and all the chapters of the 
ISSI book ??. 


That stellar variability provides deep insig hts into stellar structure has been known 
since the beginning of the 20th century (e.g. IShaDlevlll914l 'l. In nonradial oscillators 
such as the Sun, the numerous pulsation modes with distinct frequencies encode de¬ 
tailed information about their interior structure. The Sun oscillates in approximately 
10 million distinct oscillation modes and with the help of space-borne instruments we 
are now able to measure some 6000 individual modes and up to several million modes 
via so-called ridge fitting. The amplitudes of solar-like oscillations, however, are small 
compared to those observed in classical pulsators, such as Cepheids or Delta Scuti stars, 
and their frequencies can also be quite high (5 minutes in the solar case). Therefore 
only the latest generation of instruments have been able to detect solar-like oscillations 
in stars other than the Sun. 


The first indication of excess acoustic power of solar-like oscill ations in other stars 
with a frequency dependence similar to the Sun was reported bv iBrown et al.l ( 1 9911 1 
for the F5 star Procyon A (a CMi), for which th e first unambiguous detection of 
solar-l ike oscillations wa s repo rted bv iMartic et al.l 1 1999ll from the ground, as well 
as bv ISchou and Buzasil (l200ll) using data collecte d by the WIRE spacecraft. These 
results were conhrmed later bv iMosser et al.l (l2008h . The first (unconhrmed) detection 
of individual peaks in the acoustic power spectrum from high-precisio n time-resolved 
spectr oscopic observations was published for the GO star rj Boo by iKieldsen et al.l 
(1191^), but it was not until 2003 that an u nambiguous conhrmation was established 
bv JGarrier et al.ll200.ll : IKieldsen et ahllioo,^ . 


High-degree oscillation modes will not be accessible in other stars in the foreseeable 
future because they require spatially resolved variability data. However, the observable 
low-degree modes in distant stars still provide rich information about the internal 
structure and dynamics of a large number of solar-like stars at different evolutionary 
stages, as demonstrated recently by the results obtained by the CoRoT and Kepler 
missions. 


Additional p rojects, such as the Sun-in-time I Gridell [20071 1 or the Sun as a star 
1 GiamDaDall200fil ~l. try to characterize the Sun by comparing its property to solar-like 
stars at various evolutionary stages. In particular, the understanding of solar rotation, 
convection, variability and magnetic activity in a stellar context allows to better con¬ 
straint theory by deriving key scaling laws as a function of age, rotation, chemical 
abundances and mass of those various physical processes. We defer to §1.4 to 1.6 for 
discussions of the relevant observations and to §2 for the description of the current 
theoretical understanding. 


Returning to asteroseismic diagnostic, even with the limited seismic information 
of low-degree oscillation modes, some theoretically important properties have become 
accessible, such as the gross structure of the energy-generating core and the extent to 
which it is convective, and a large-scale core to envelope contrast i n the angular velocity 
of evolved stars (e.g. iMosser et al.l I 2 OI 2 I : IPeheuvels et al.l 120121 : iMathur et al.l 120121 : 
iGarcfallioisl . and references therein). Such information will be of crucial importance 
for checking, and then calibrating, the theory of the structure and evolution of stars. 
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One of the many interesting applications of utilizing low-degree modes for mea¬ 
suring various aspects of the structure in distance stars is the measurement of rapid 
variations in the stellar background, brought about by, for example, the rapid variation 
in the acoustic cutoff frequency at the base of the surface convection zone or the vari¬ 
ation in the first adiabatic exponent in the helium ionization zones. The design and 
analysis of low-degree seismic signatures of the structure of spherical stars is there¬ 
fore of central importance to asteroseismic diagnosis, and va rious research groups have 
been working on improved seismic diagnostic techniques fe.g. lHoudek and Goughll2007l . 
and references therein). First applications of low-degree seismic diagnostic techniques 
have been carried out r ecently on 19 solar-type stars observed by the Kepler mission 
( Mazumdar et ^|2014| ~). 

However our poor understanding of the physics of the stellar surface layers causes 
a major problem in the analysis of global solar-like oscillation, because the proper¬ 
ties of the near-surface layers have a major impact on oscillation frequencies. Typ¬ 
ical stellar structure calculations treat the superficial layers with simplified atmo¬ 
sphere models and the stratification of the supera diabatic region by m eans of a lo¬ 
cal, time-independent mixing-length approach (e.g. iBohm- Vitensel Il958l l . The effects 
of the (turbulent) Reynolds stresses in both the hydrostatic equilibrium and pulsa¬ 
tion calculations are essentially always ignored. Linear pulsation calculations typi¬ 
cally adopt the adiabatic approximation and ignore the momentum flux perturba¬ 
tions (turbulent pressure fluctuations). However, nonadiabatic effects and the fluctua¬ 
tions of the turbulent fluxes (heat and momentum) do modify the mod elled pulsation 
eigenfunctions and consequen t ly als o the oscillation frequencies (e.g. iGoughI Il984al : 


iBalmfoi^^ 199^ Houde^ 1996 j_ 201^ . For a “standard” solar model, such as Model S 
IchristensemPalsgaarT'^*^ 199fil i. the frequency differences, scaled with the mode 
inertia, between the solar and model frequencies are dominated by the near-surface 
effects and are predominantly a function of frequency alone. Differences between solar 
frequencies observed with the GONG instruments and adiabatically computed Model 
S frequencies are illustrated by the the symbols (plusses) in the left panel of Figure 1 
(colours indicate different spherical degrees). The increase of the frequency residuals 
with oscillation frequency depends on the modelling details of the functional form of 
the acoustic cutoff frequency (as it affects the acoustic potential) with radius in the 
near-surface layers. In the Sun the acoustic cutoff frequency is about 5.5 mHz. 

Sophisticated 3-D hydrodynamical simulations of stellar convection have been used 
to estimate the effect of turbulent pres s ure in the equilibrium model (the mean Reynolds 
stress). For example. iR.osenthal et all ( 1995h investigated the effect on adiabatic eigen- 
frequencies of the contribution that the turbulent pressure makes to the mean hydro¬ 
static stratification. They examined a hydrodynamical simulation bv istein and Nordlundl 
(Il99ll') of the outer 2% by radius of the Sun, matched continuously in sound speed to a 
model envelope calculated, as in a “standard” solar model, with a local mixing-length 
formulation. The resulting frequency shifts of adiabatic oscillations between the simu¬ 
lations and the “standard” solar reference model, Model S, are illustrated in the right 
panel of Figure 1. The frequency residuals behave similarly to the data (plus symbols) 
in the left panel of Figure 1 though with larger shifts at higher oscillation frequencies. 

The effects of nonadiabatic ity and convection dynamics on the pulsati on fre quencies 
were, for example, studied bv iBalmforthl (Il992l l. fRosenthal et al.l lll995ll and iHoudekI 


(119961 1. In these studi es the n onlocal, time -dependent generalization of the mixing- 
length formulation bv IGou^ ( 1977alfl977bh was adopted to model the heat and mo¬ 
mentum flux consistently in both the equilibrium envelope model and in the linear 













































Fig. 1 Figure 1: Frequency residuals between the Sun and solar models as a function of 
frequency. Left: scaled differences between observed GONG frequencies and adiabatically com¬ 
puted frequencies of the “standard” sola r Model S. Different colours indica te different degrees 
of the oscillations modes (adapted from IChristensen-Dalsgaard et '^1199611 . The curves show 
frequency residuals between models computed with a nonlocal time-dependent and a “stan¬ 
dard” mixing-length formulation for convection: The solid (blue) curve shows the adiabatic 
residuals cause by the Reynolds stress in the mean model. The dashed (red) curve is the fre¬ 
quency shift caused by nonadiabaticity and conve ction dynamic s. The overall frequency shift 
is plotted by the dot-dashed curve (adapted from [Houdekl[T996h . Right: Scaled adiabatic fre¬ 
quency differences between a model for which the near-surface layers were represented by a 
hydro dynamical simulation, and the “standard” solar Model S (adapted from I Rosenthal et al.l 
Il995f) . 


nonadiabatic stabi lity analysis. The ou tcome of such calculations is shown in the left 
panel of Figure 1 ( Houdekl ll99fil . l20icll by the dot-dashed curve. Similar as for the 
hydrodynamical simulations the effect of the Reynolds stresses in the mean structure 
decreases the adiabatic frequencies (solid curve) with a maximum deficit of about 
12/rHz, which is much smaller than the deficit of about 24^Hz from the hydrodynam¬ 
ical simulations (symbols in the right panel). 

The effects of nonadiabaticity and convection dynamics (dashed curve in the left panel), 
however, lead to an increase of the mode frequencies by as much as ~9^Hz, nearly can¬ 
celling the downshifts from the effect of the Reynolds stress (turbulent pressure) in the 
mean structure (solid curve in the left panel), leading to the overall effect indicated by 
the dot-dashed curve in the left panel. If the positive frequency shifts between mod¬ 
els computed with and without nonadiabatic effects and convection dynamics (dashed 
curve in the left panel) are interpreted as the nonadiabatic and momentum flux correc¬ 
tions to the oscillation frequencies then their effects are to bring the frequency residuals 
of the hydrodynamical simulations (right panel) in better agreement with the data (plus 
symbols in the left panel of Figure 1). If we were to understand the reason why the 
semi-analytical, nonlocal convection model underestimates the effect of the Reynolds 
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stress in the mean model on the adiabatic eigenfrequencies, we would come closer to a 
self-consistent explanation for the still unsolved problem of reproducing theoretically 
the observed high-frequency oscillations in the Sun and in sun-like stars. Currently only 
simp le procedures for estimating the near-s urface frequency cor rections ( Kieldsen et ahl 
I 2 OO 8 I I are adopted by the community fe.g. iMathur et al]l2012h . This simple procedure 
is based on the calculati ons of a model, in which the frequency s hifts are scaled by a 
power-law in frequency ( Christensen-Dalsgaard and Goughlll980l ). A physically more 
realistic description of the surface effects with reliable applications to solar-like stars is 
therefore of imminent importance for the whole community to determine more reliable 
radii and ages in distant stars from the observed mean large frequency separation and 
small frequency separation, for they depend crucially on the near-sur face frequency 
corrections fe.g. IChristensen-Dalsgaard and GoughiriQSOl : lHoudekll201oll . 


1.3 Asteroseismic Scaling Relations 


With the advent of space missions, asteroseismologists studying solar-like stars started 
to apply techniques originally developed for Sun-as-a-star observations to other stars. 
With only one or two solar-li ke stars per target field observed by GoRoT (e . g. HD 49933, 
HD 181420 and HD 181906. [Appourchaux et al.ll2008l : lBarban et al.ll2009l : lGarcfa et all 
l2009l . respectively), it was possible to do a full analysis of the power spectrum, deter¬ 
mine the individual p-mode freq uencies and, after combining the seismic observables 
with the spectroscopic ones fe.g. iBrunttI 2009ll. perform a full modeling of each target 
using stellar evolution codes (e.g. Mathur et al.l[^13bh . Later on, with the study of 
red giants in the GoRoT exoplanet field and the arrival of Kepler, asteroseismologists 
faced a new challenge: the need to analyze massive samples of thousands of stars with 
high-quality light curves. The combination of the high data volume and the surface 
effects above made manual analysis impractical. As a result, seismologists turned to 
characterizing the pulsation frequency patterns with a measurement of the typical fre¬ 
quency of oscillation and the typical frequency spacing between carefully selected and 
classified modes. To do this for large sample, automatic pipelines were developed in 
order to extract these mean seismic observables: Av and Vmax (e.g. Huber et al.ll200^ : 


iMosser and AppourchauJl2009l : [Hekker et "^120101 : iMathur et al.ll^l0a l. 

The large frequency spacing, Aviin) = i is proportional to the sound 

travel time in the cavity in which the mode propagates and, therefore, it is proportional 
to the density inside the star : Av ~ (M /. This relation shows a very good 
agreement with modeling (e. g. lUlric hl ll986h and can be scaled from the solar values as 
( Kieldsen and Beddinel[l99!Tl l: 


Av « Av(t, 


1/2 


-3/2 


( 1 ) 


M \ ' ( R 
Mq) [Rq 

in which ZiiZQ = 135.1 ± 0.1 /iHz, as derived bv iHuber et al.l 1 201 il l using 111 VIRGO 
subseries of 30-day, each spanning fro m 1996 to 2005 and analyzed in the same way 
as the Kepler data. More recent work ( White et al.ll2oi3 : lMosser et al.ll2013l l based on 
numerical si mulations sug gests the need for small correction terms. 

In 1991, iBrown et al.l suggested that the frequency at the maximum power of 
the p-mode bump, z zmax, should sc ale with the acoustic cutoff frequency (see also 
iBelkacem et al.l 120111 1 , and following iKieldsen and Beddind ( 1995h it could be scaled 
from the solar values according to: 
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^max 


i^max,0 


M \ 
M^J 



( 2 ) 


in which = 5770 K, and i-'max,© = 3090 ± 30 /iHz I Huber et al ]H). This 

constraint is less strongly motivated by theory and requires empirical guidance. 

Therefore, from equations m and m, one can obtain the stellar mass and radius 
as a function of the seismic variables i/m ax and Av - assum ing that Tgg is known 
independently of evolutionary models (e.g. iMiglio et al.ll2009ll : 



_ It has also recently be en shown le.g. lStello et al.l2009l : [Hekker et "SlbooglJMosser et al.l 

l2O10l : lHekker et al. IS) that solar-like oscillations in main-sequence stars follow to a 
good approximation the following relation: 


Av 




t^max 

^max,© 


( 5 ) 


where b has often values around 0.75. Figure[5](top) shows the relation Av versus Vniax 
for 1700 stars from the main sequence phase (black diamonds at the right hand side) to 
the red clump phase (red triangles appearing in a diagonal branch between 20 to 50 pHz 
in the bottom panel) observed by the Kepler satellit e. Although the relation appears 
to be constant for all of these stars, several authors ( Mosser et al.l 12010 : iHuber et al.l 
l20inll have suggested that the slope is different between red-giant and main-sequence 
stars. To enhance such a difference, we have subtracted the luminosity dependence by 
raising r/max to the power of b = 0.75. A fit to the residuals below and above t'max= 
300 /iHz - which roughly marks the transition from low-luminosity red giants to sub¬ 
giants results in a steeper slope with increasing i^max. It is important to note that for 
i^max close to the solar value, for example, the use of a power-law relation calibrated 
to red-giant stars would lead to an underestimation in Av by ~ 10%. 

Very precise surface-gravity values, {g), - with uncertainties much less than 0.05 
dex - can also be inferred directly from the global seismic parameters and from the 
effective temperature: q « (e.g. iMorel and MigliJ I 2 OI 2 I : iThvgesen et al.l 

12012 ! : lilekker et al.ll2013h . This is very important because the seismic logg values can 
be used as an independent fixed parameter i n the spectroscopic analy sis to compute 
the effective temperatures in an iterative way i Silva Aguirre et al.ll20ial . Indeed, both 
parameters are correlated and cannot be determined independently from the spectra. 
Therefore, the combination of a seismic loggr and a spectroscopic Tgg can reduce the 
uncertainties in both parameters. CoRoT and Kepler stars with asteroseismic surfac e 
gravities will be used as calibration stars for the GAIA mission ( Creevev et al ][2^. 
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Fig. 2 Av versus i^max (top panel) for 1700 stars observed by Kepler taken from iHuber et all 
1 I 2 OIII I . The red triangles show stars observed in long cadence (red giants), while black diamonds 
are stars observed in short cadence (solar-like stars). In the bottom panel we have shown the 
same relation but after removing the luminosity dependence by raising i^max to the power of 
0.75. Green lines show power law fits to the Au — i^max relation for two different intervals 
of fmax- The blue d ashed line shows t he relation of Eq. [S] derived using both red giants and 
main-sequence stars JStello et al.ll200^ 1. 


1.4 Stellar Convection and the Background Signal in Oscillation Spectra 


Solar-like pulsating stars are those in which their acoustic mo des are excited by tur¬ 
bulent motions occurring in their convective outer layers (e.g , 

I 1977 I : lBalmforthlll992l : iHoudek et al.l[l99^ : IChanlin et al.ll2005 


Goldreic^^ndKeeleyl 

Beikacei^^r*^ L 


2008 


The first observations of granulation are quit e old jHerschell 1801 ). but it is not until 


1933 that their turbulent origin was unveiled I SiedentoDifT93^ )~ 

In solar-like stars, exhibiting outer convective envelopes, the granules evolve with 
time producing small photometric fluctuations. Therefore, the analysis of stellar light 
curves in the Fourier domain shows a continuous “background” on top of which the 
acoustic modes are superimposed (see some examples of main sequence stars, subgiants 


[Ballot et al.ll20l'il : iMEdhureT*^ 201 Ibl : Metcalfe’^Tlin 20121 : iMathur^TliL 2013bl ). 


and red giants in Kallinger et al.l 2O10l: Chaplin et a .1 201(1: Camnante et al.l 201 ll : 


The first successful attempt to model this continuous background was done by 
iHarvevI ( 1985h . He approximated each scale of the convective motions (granulation, 
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supergranulation...) by an exponentially decaying time function corresponding to the 
autocovariance of its time evolution with two parameters: the characteris tic time of the 
granul es and the total brightness fluctuations associated with them Isee lMathur et ahl 
l2011al . for a complete review of the different formulations commonly used in astero- 
seismology). Moreover, it has been demonstrated very recently that the properties of 
the granulation and the global seismic variables are related (e.g. IChaolin et al.ll2011bl : 
iMathur et ahllioilah . It is also important to notice that in some stars another feature 
is visible close to the p-mode envelo pe that has been related with the presence of fac- 
ula e in the star (Karoff et al.| 1^131 . and references therein). Recent theoretical work 
by ISamadi et al ] (l201'ji has also tried to explain the observed relationship between 
granulation amplitude in the oscillation spectrum and stellar spectral type. 

In the context of stellar dynamics, putting observational constraints in the convec¬ 
tion inside stars is very important because convection is one of the major ingredients to 
understand stellar magnetism and dynamo processes. Therefore, precise time domain 
measurements can offer direct constraints on the granulation characteristics as well as 
the precise determination of the convective outer envelope. Recently, using Kepler data, 
it was for example possible to calibrate the mixing-length parameter omlt used in clas¬ 
sical 1-D treatment of convection , and to show that t he solar value, us ually used for 
gener ating grids of stellar models ( Ekstrom et al ]H) is overe stimated ([B qnaca et al.1 
I 2 OI 2 II . Likewise granulatio n data can put constraints on local dMaeic et al.ll201^ and 
global ( Miesch et al. llili stellar convection simulations, that could in turn help us to 
characterize stellar oscillation s, and convection parameters such as the mixing length 
I Tramnedach and Steiniboilh . This important by-product of asteroseismology studies 
can help to better constrain convective heat transport in 1-D stellar models, and hence 
to constrain the stellar radius which depends on omlt- 

Asteroseismology can also provide important information on the depth of con¬ 
vective envelopes. At the base of the surface convection zone, there is an abrupt 
change of stratification (going from adiabatic to subadiabatic) which impacts the 
sound speed and density profiles (see the pionnerin g work on the base of the so¬ 
lar convection by I Christensen-Dalsgaard et al ] ll99lli . By using appropriate seismic 
diagnostic techniq ues it is then possible to infer the extent of the convective zone 
of so lar-like stars I Goughiri99Cll : iBallot et al]l2004l : lHoudeklEo04l : iHoudek and GoughI 
[ 2007 ). This approach h as now been successfully applied to Kepler data for field stars 
dMazumdar et^l2014l i. 

Complementary to the analysis of asteroseismic data, numerical simulations to 
model and describe the turbulent state of stellar convective enveloppes or cores have 
been pursued. Several groups have undergone a systematic study of the global prop¬ 
erties of stellar convections and their associated mean flows (differential rotation and 
meridional circulation). We refer to §2 for a more detailed discussion of numerical 
simulations of rotating and magnetized convective stars. 


1.5 Stellar Rotation and Spin Down 


All stars rotate, and understanding their rotation history is key to our overall under¬ 
standing of their dynamics. The angular momentum evolution of stars is an initial value 
problem, and a wide range of rota tion rates are observed in the youngest protostars 


(see for example iRebull et al.l 20061. for the Orion Nebula Cluster). The relatively slow 


rotation of such stars ( Vogel and Kuhilll98ll ') demands efficient angular momentum loss 
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from the infalling material du ring the hydrod ynamic collapse phase, probably through 
strong magnetic interactions ( Slni_et_an 19^) . Gaseous accretion disks persist around 
protostars for up to 12 Myr ( Bell et al. IM), and the exchange of angular momentum 
between them and their disks can dramatically modify their evolution. In the limit of 
strong magnetic interaction, star-disk coupling can regulate stellar rotation, preven ting 
protostar s from spinning u p as they contract ( Koenig]lll99ll : iKepoens et al.fll99Sh but 
see also ( Matt et al. IS). The early angular momentum evolution of stars is there¬ 
fore a window into the star and planet formation process. Stars also reach the main 
sequence with a range of rotation rates. T his can cause star-t o-star differences in the 
degree of rotationally induced mixing (see IPinsonneaultlll997l . for a review). Once on 
the main sequence, the surface rotation of stars responds both to internal angular mo¬ 
mentum transport and to angular momentum loss from magnetized solar-like winds. 
There are two strong practical applications that result: 

— 1) We can use the observed angular momentum evolution as a test of stellar physics. 
Many instabilities and physical processes such as dynamo action, mixing, large scale 
flows are directly connected to the rotational rate of stars. Stars differ in their ther¬ 
mal structure, rotation rates, and lifetimes; we therefore have natural laboratories 
for testing physical models. This is especially important for our understanding of 
solar-like winds, as angular momentum loss driven by them is a central driver of 
angular momentum evolution in low mass stars. 

— 2) With calibrated models we ca n then use rotation as a stellar population di¬ 
agnostic. IWeber and Davi3 ( 1967h explained the slow solar rotation as a conse¬ 
quence of angular momentum loss from a magnetized solar-like wind . They pre¬ 
dicted dJ/dt ~ which implies f2 ~ in the asymptotic limit. ISkumanichl 


( I 972 II confirmed this time dependence in rotation, chromospheric activity, and 
lithium depletion in open cluster stars. Rotation is therefore a potential age in¬ 
dicator, provided both that it can be measured and that the rotation-mass-age 
relationship can be calibrated. 

Both our observational knowledge and the sophistication of our theoretical models 
have advanced considerably since this pioneering work. We therefore briefly summarize 
the observational picture (see Figure 3) before proceeding to the theoretical picture 
that has emerged. Rotation rates can be inferred spectroscopically from Doppler line 
broadening. Such measurements, however, suffer from projection effects (usini is what 
is measured) and become difficult when the rotational broadening becomes compara¬ 
ble to or smaller than other sources, such as microturbulence or thermal broadening. 
Rotation velocity measurements are therefore primarily used for young and massive 
stars, and through the 1980s and early 1990s most ro tation data was spectroscopic. 
In a pioneering paper, Istauffer and HartmannI (Il987h demonstrated that cool stars 
in the Pleiades arrived on the main sequence with a wide range of rotation rates, 
which created three significant problems for simple angular momen tum evolution mod¬ 
els. E xtremely large torques would be predicted for rapid rotators ( Pinsonneault et al.l 
Il99fll 'l. and an unsaturated wind law would not permit the survival of high rotation 
rates onto the main sequence. Magnetized solar-like winds are too inefficient to spin 
down young stars to the observed rates for the slowest rotators, requiring an addi¬ 
tional mechanism for extracting angular momentum in the pre-main sequence. Finally, 
once on the main sequence, there is a rapid initial drop in rotation rate, followed 
by a pause where surface rotation changes little; older stars then resume spinning 
down at the predicted asymptotic rate. This is evidence for core-envelope decoupling 
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( MacGregor and BremieiJ ^ 199l];_ Krjghj 2 ^murthj_et^^ 1997 ) on timesc ales of tens to 
hundreds of Mvr fsee lDenissenkov^niLll20innGaUerand Boxivierll201,fl . for recent dis¬ 
cussions). The observed pattern is therefore richer than a simple power-law relationship, 
requiring both more sophisticated theory and stronger empirical constraints. 

Star spots cause periodic brightness modulations as they transit across the visible 
disk. Dedicated time domain studies of stars can therefore measure rotation periods. 
Rotation period measurements do not suffer from inclination effects (although they 
require a favorable inclination angle, e.g. not pole-on, to produce a signal.) There are, 
however, significant observational selection effects that favor the detection of short pe¬ 
riod systems, which also tend to be high amplitude. There has been an enormous 
increase in the quantity and quality of stellar rotation period data. In particular, 
the advent of large surveys designed to search for extrasolar planets has yielded very 
large data sets of ste llar rotation periods as a natural byproduct for both field stars 
dHartman et al.l[201ll ) and open clusters stars . Exa mples of the latter in clude the exten¬ 
sive MONITOR program ( Aigrain et ahllioOTl ) fsee lMoraux et al1l2013l. for the mo st re¬ 
cent results), t he SuperWASP data f o r the Hyades and Praesepe ( Delorme et al ][m3) 
and the large I Hartman et al ] l2009l . l20inl) M37 and Pleiades data sets. The Kepler 
and CoRoT missions reach lower amplitudes and can detect stars with longer rota¬ 
tion periods, pushing ste l lar rotation measurements into the solar regime dSasri et al.l 
2011 ; Garcia et al. 201^ Lanza 20^ Walkowicz and Basril 1201.31 : iNielsen et al. 201.31 : 


Reinhold et al .Il201.3l^ McQuillan et al ]|20 i 4) . With the aid of these samples, a standard 


empirical picture of angular momentum evolution has emerged, with the following basic 
ingredients. 

Stars are born with a range of rotation rates, and they would be expected to spin 
up as they contract and approach the main sequence. Through protostar-disk inter¬ 
actions some stars lose additional angular momentum in the first ~10 Myr, widening 
the observed distribution. The slow rotator population cannot be generated from the 
young star for ming regions, as not ed above, but is already in place in the 13 Myr old 
system H Per ( Moraux et al. iim. This confirms the hypothesis that protostar-disk 
interactions are important for stellar angular momentum evolution. The distribution of 
rotation rates is initially broad and narrows with age, which is a natural consequence of 
a torque that increases more than linearly with rotation rate. This convergence process 
is most rapid for solar mass stars and the timescale for spin down increases with de¬ 
creased mass. This is usually interpreted as evidence for mass-dependent saturation in 
the angular momentum loss rate, typically parameterized with a Rossby scaling (rota¬ 
tion period relative to the convective overturn timescale.) Core-envelope decoupling is 
requir ed to explain young slow rotators, with a time scale that increases with decreased 
mass ( Denissenkov et al.lf20ld : lDelorme et al Jim. There is no break or discontinuity 
in the observed spin down behavior at the fully convective boundary, which is surpris¬ 
ing given dynamo model predictions. New data sets in older solar analogs in Kepler 
stars will provide tests of the mass and composition dependence of rotation in old 
dwarfs, which is currently lacking. Field star studies will need to consider the presence 
of background populations when modeling angular momentum evolution. For example, 
angular momentum lost from tidally synchronized binaries is extracted from the or¬ 
bital rather than the spin angular momentu m reservoir, which resu lts in very different 
angular momentum evolution properties f see I Andronov et al.ll200^ . for a discussion in 
the context of blue stragglers). 

These large data sets have spurred rene wed interest in the usage of stellar rota¬ 
tion as an age indicator, or gyrochronology [BarnesI ( 20031 ). Distinct theoretical seal- 
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Fig. 3 The observed distribution of rotation rates in stars of dif f erent ages. Data for the 
young protostars in H Per (upper left) is taken from iM^oraux et al.l ll2013lb The data for the 
125 Myr old Pleiades system (lower left) is t aken from the reanalys is by Coker Sz Pinsonneault 
(2 015, in prep). M37 data is orig i nally from iHartman et al.l ll2Q09ri as interpolat ed onto mass 
by Ivan Saders and Pinsonneault! <l2013lb Data for old field dwarfs is taken from [Garcia et all 

il2ni4h . 


itiRS and fittinR function have been proposed (e.R. Mamaiek and Hillenbraridl l2008l : 

' 312 


iBarnes and KinJl2010l : lEpstein and Pinsonneauitll2C)14 , for a discussion). There are sig¬ 

nificant zero-point shifts between these formulations, which will require resolution when 
compared with empirical data. There are strict lower bounds t o the uncertainties in gy- 
rochro nology relations even for an ideal theoretical framework. lEpstein and Pinsonneault! 
(l2014h identified two fundamental limiting effects for the precision of gyrochronology: 
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The intrinsic range of rotation rates at fixed age and the range of surface rotation 
rates in a given star from latitudinal differential rotation. The former is more im¬ 
portant for young and low mass stars, while the latter dominates for older higher 
mass stars. The net effect i s an a ge uncertainty at the 2 Gyr level for old stars. 
Ivan Saders and Pinsonneaultl 1 20131 1 extended this to subgiants, and argue for a differ¬ 
ent gyro chronology relationship for them than for dwarfs. Recently the gyrochronology 
concept has been adapted to consider magnetic f ield instead of rotation as a proxy for 
age: e.g. magnetochronology dVidotto et al.ll2014h . It is indeed clear (see next section) 
that stellar magnetic properties vary over t he course of stellar evolution, older s tars 
being less active their younger counterparts ( GiamDapalboOiTl : lOregorv et al.ll2ni2ll . 

The theoretical interpretation, numerical models and multi-D simulations allowing 
us to study and to understand the stellar age-rotation-activity relationship will be dis¬ 
cussed in S2. 


1.6 Stellar Activity & Magnetism 


The great advantage of comparing the Sun to other stars is to disentangle what is 
specific to our star from what is generic to solar-like stars. Systematic observations of 
rotation and magnetic properties of solar-like stars have revealed interesting trends. 
For instance surface differential rotation is fo und to increase with Tpff; F-star possess 
a larger latitudinal con trast than K-stars do ( Donahue et al.|[l996l : iBarnes et al.ll2005l : 
ICollier Cameronll200it). Similarly a rotation- activity relationship has been known for 
quite some time now ( Pizzolat o et al.l 200311. A relation between cycle and rotation 
periods has also been advocated ( Noves et al.lll984bl : iBohm-Viten^ 2007ll . 

Meaningful information regarding magnetic cycles underlying stellar activity, or 
alternatively to rule out the existence of magnetic cycles, require long-term stellar 
activity observations and their analysis. Direct observations of starspots is a relatively 
new capability and studies of their behavior over timescales long enough to be relevant 
for stellar dynamos, do not yet exist. Nonetheless, indirect evidence of magnetic activity 
cycles (or the lack thereof!) can be obtained by different techniques. Our traditional 
tools for measuring stellar cycles have been synoptic stellar chromospheric activity 
observations, complemented with more limited stellar coronal X-ray variability data. 
Asterseismology may prove a useful additional set of data. 

The most important contribution in the context of long-term stellar activity obser¬ 
vations has been the seminal Mount Wilson Observatory (MWO) Ca lcium (Ca) H-l-K 
Project ( Wilsoniri978bl : lNoves et ai1ll984bl : lBaliunas et al.lll98!Tl . ll995|i . The emission in 
these line cores, expressed as a Ca H-fK activity index, is understood to be a measure 
of the magnetically mediated, non-thermal heating o f the chromosphere (provided that 
the b asal thermal flux is properly accounted for, e.g. ISchriiver et al.iri989f : ICuntz et al.l 
I 1999 I I. Thus, the long-term variation of the Ca H-l-K in dex can be taken as a proxy for 
variability in stellar magnetic fields (see the review bv lHal]ll2008ll . The Lowell Obser¬ 
vatory’s Solar- Stellar Spectrograph program provides comp lementary observations of 
Sun-like stars I Hall et ^ I 2 OO 7 I : lHallll20(Sl : iHall et ^l2009l l adding to the knowledge 


base of the MWO program. An independent measure of stellar magnetic output is 
the coronal X-ray flux - which originates in the magnetic heat ing of stellar corona e 
and which is even more tightly correlated with magnetic flux ( Pevtsov et al. IH). 
Although stellar X-ray cycle observations are rare, the amplitude of stellar X-ray flux 
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have been observ ed for multiple stars i Hemoelmann et al.lll99^ : iMicela and Marinol 
I 2 OO. 1 I : [Oudellbooll and this usefully constrains the strength of the dynamo generated 
magnetic fi elds. More recent su rveys analyzing the a c tivity of solar-like stars are to 
be found in lWright et al.l ( 2004l ~) and iGiampapa et ahl (l2006l i. We provide here a brief 
summary of the important trends and relationships relevant for stellar activity cycles 
that have been gleaned from these long-term observations. 


1.6.1 General Activity Trends with Rotation 


First we focus on the scaling behavior of magnetic activity with the stellar rotation 
rate. Stars that are faster rotators have higher magnetic activity amplitudes than 
slower rotators. They are usually variable, but do not display regular cycles. Stars 
with intermediate rotation rates have intermediate magnetic activity output, which 
sometimes display periodic behavior. Stars which rotate relatively slowly, such as the 
Sun, are more likely to have regular, periodic cycles, but their activity amplitudes 
are also relatively smaller. The latter category also sometimes display “flat” (i.e., non¬ 
varying) activity with low amplitude, akin to global magnetic minima such as the solar 
Maunder minimum. We note here that fast rotators also sometimes display a “flat” 
(non-varying) activity profile albeit with high activity amplitudes- which is not to be 
confused with global minimum like states. 

In general the activity amplitudes increase with rotation r ates, but for very fast 
rotators, the measured amplitudes tend to reach a saturation I Pizzolato et al.ll2003l : 


IWright et al.ll201lh . which in X-rays is of order \og{Lx /~ —3.2. For G type stars 
this saturation is found for rotation rate above 35 kms“^, for K type stars at about 
10 kms“^ and for M dwarfs around 3-4 kms“^. How stellar magnetic flux scales with 
rotation rate is important to understand, since it is telling us how the magnetic field 
generated by dynamo action inside the stars emerges and imprints the stellar surface. 
Understanding the apparent saturation phenomenon is also i mportant : for e xample, 
it could arise if only the filling factor / reaches a maximum ( ReinersI l201j ). Multi- 
perio dic cycles are also sometimes observed, especially in faster rotators 1 Olah et al.1 
l2009l i. For those stars with cyclic magnetic activity, the systematic analysis of stellar 
data revealed that for solar type stars there is a good correlation between the cycle 
and rotation periods of the stars and that correlation is even stronger when using the 
Rossby number (Ro = Protlr) that takes in to account the convection turnover time r at 
the base of the stella r convective envelope i Noves et al.ll 198^11984^ : ISoon et aklllfi^ : 
iBaliunas et al.iri99fih . As the star rotates faster, its cy cle period is found to be shorter. 
This trend is summarized in Figured More precisely, NoYgs_etL_al lll984ah found that 


Pcyc oc ProtJ with n = 1.25 ± 0.5. Saar and Brandenburg 1 199£ 1: Isaail ( 2002l i have 


used an extended stellar sample to argue that there are actually two branches when 
plotting the cycle period vs the rotation period of the stars: the primary (starspot) 
cycle and the Gleissberg, or grand minima type, modulation of stellar activity. For the 
active branch they found an exponent n ~ 0.8 and for the inactive one they inferred 
n ~ 1.15. 


1.6.2 Maunder Minima Stars: 

Many of the stars (roughly 15%, see lGiampapa et al.ll200t^ l. studied in the Mount Wil¬ 
son surveys, displayed activity that were flat (i.e., non-cyclic) and at or below current 
solar minimum levels. Similar studies with Kepler data recently found that among 
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Rossby Number 



Fig. 4 Stellar activity trends as determined from observations of Ca H+K and X-ray emission 
from F-G-K stars. The arrows in the axis label show the direction of increase of the rotation 
rate and Rossby number. In general, magnetic cycle amplitude increases while cycle period 
decreases with increasing rotation rate (i.e., decreasing rotation period). The critical reader 
is referred to the original o bservations cit ed in the text for the exact scaling laws and other 
subtleties in the data. From lN^ d3 <l2n04l'l . 


the sampled solar-like seismic stars about 2% show a flat activity level. One possi¬ 
ble interpretation of these results is that these are “Maunder minima” stars where 
the global magnetic activity has temporarily ceased. Complement ary studies of coro¬ 
nal emission also point to th e existence of such low-activity stars dPoppenhager et al.l 
I 2 OO 9 I I . iBaliunas' et al.l dlQQSh also provides evidence indicative of a star entering into a 
Maunder minimum like phase - with its activity variation flattening out accompanied 
by a reduction in amplitude. 


1.6.3 Starspots, Polar Caps and Stellar Butterfly Diagrams: 

Recent advances in stellar photometry and doppler imaging techniq ues now allow us to 
observe starspots and their dynamics. Such observations (reviewed in lBerdyuginalliooll 
show that there is oftentimes a preferred longitude for starspots. Faster rotators tend 
to have large polar spots which contributes to a strong polar cap, while in general, in 
very active stars, the spot distribution tends to be more uniform across all latitudes. In 
some stars, the spots evidently appear at mid-low latitudes and subsequently migrate to 
higher latitudes, reminiscent of the sunspot migration (although towards low latitude) 
that generates the solar butterfly diagram. These observations do not span a long period 
of time and hence detailed studies of stellar butterfly diagrams is still in its infancy. 
Nevertheless, observations of spot emergence latitudes and migration are expected to 
provide more stringent constraints in the future for stellar dynamo and flux emergence 
models. 
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1.6.4 A stem seismology and Activity 


Ensemble asteroseismology can now add to our developing understanding of stellar ac¬ 
tivity obtained by other methods, such as spectropolarimetry and Call H & K studies. 

The first contribution of helio and astro-seismology comes through their ability 
to provide stellar properties and, in particular, stellar ages for field stars that are 
otherwise difficult to analyze. Our lack of knowledge of stellar ages has been an im¬ 
portant limitation in our understanding of the Mt. Wilson survey data, for example. 
Many of the solar-like stars showing a fiat S-index in the original sample analyzed by 
iBaliunas and .Tastrowl dlQQfll l -and assumed to be in a similar state to the Maunder 
minimum in the Sun- were found to be more evolved stars (which would be expected 
to less a ctive) when a ccurate parallaxes from the Hipparcos mission were taken into 
account (IWrightll2004h . Therefore, accurate stellar properties are essential for proper 
int erpretation of their magne tic indexes (see an extended discussion in the review 
bv l.ludge and ThomDsonll2012h . But seismology can also be employed to directly study 
stellar activity. Magnetic fields modify the conditions inside the stars and the properties 
of the modes are changed accordingly. Seismology can therefore provide complemen¬ 
tary information to the one obtained by the above mentioned techniques. Moreover, it 
can for example provide evidences of the existen ce of subsurface on-going magnetic cy¬ 
cles even during Maunder minimum-like states ( Beer et al.ll 19981 ) in which no external 
manifestation of the internal magnetism can be observed. Howe ver, the amplitude of 
the m odes decreases for active stars as observed in the S un (e.g. I.Timenez-Reves et al.l 
I 2 OO 3 I ). HD 49933 observed by Co RoT ( Garci'a et al.l[201ol ). and solar-like stars observed 
by Kepler ( Chaplin et al.ll2011al ) difficulting their observation. 

The first diagnostic or “proxy for magnetic activity” used in helioseismology was 
the frequency shift of the acoustic modes, i.e., the displacement of the modes towards 
higher frequencies as the cycle evolves. These frequenc y shifts -which follow an in¬ 
verse mode-mass scaling ( Libbrecht and Woodardlll9^ )- are explained to arise from 
changes in the outer layers of the Sun along the 11-year magnetic cycle (Schwab e cycle ) 
(IColdreich et al.l 1 199 if : iBalmforth et al.l 1 19961 : IPziembowski and Goodel l2004 120051 ) . 
Therefore, the high-frequency p modes - haying higher upper turning points than 
the low-frequency p modes - are then more sensitiye to the perturbations induced by 
the magnetic field in the outer part of the Sun. Thus, these high-frequency p modes are 
obseryed to haye large frequency shifts with the Schw abe cycle, while the low-frequency 
p modes present little or no yariation with it (e.g., 


Il998l : iGelly et'al]l2002l : iHowe et al.l l2008|: 


■Tain et al 


been found in the F star HD 49933 ( SakherteT'^ 2011a ) . 


Ronan et al. I 1994 I : iGhaplin et al.l 


200 s). A similar behayior has also 


I 2 OO 8 I: iBroomhall et a 

.1 2012 !). all the rest of the solar low-deg 

ree p-modes properties 

change with time (e.g. 

Chaplin et al.ll2000l:IJimenez-Reves et al. 

l2n03l.l2007l:ISalabert et al.l 


Schwabe cycle. Another shorter timescale has been recently confirmed using helioseis¬ 
mology. Indeed, the quasi-pe riodic modulations of 1 to 2 years obseryed during the last 
twenty years in the Sun (e.g. Benevolenskaya|[l995l : iMursula et al.l[2003l ) has then been 


measured ( Fletcher et al. \mm in the frequency shift s of low-degree mod es measured 


by the Global Oscillation at Low Frequencies (GOLF [G.ffiriel_et_aL 1995|}}Jnstrument 
on board the Solar and Heliospheric Obseryatory (SoHO Do^^^^^^^3[l995]}jas well 


with the Birmingham Solar Oscillation Network ('BiSON lGhaplinTtaL[ll996l ). Al¬ 
though the origin of this modulation is not well understood yet, the latest results priy- 
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iledge the hypothesis that th e origin of the perturbation is located in the subsurface 
layers ( Simoniello et al1l2013h . 

The continuous photometric monitoring of hundreds of thousands of stars by CoRoT 
and Kepler provide long baseline stellar variability samples of unprecedented scale and 
precision. Their duration is also substantial, comprising a period of more than 4 years 
in the case of Kepler; this is long enough to begin to see the effec t of activity cycles . The 
first ensemble ana lysis Kepler light curve varibility was done bv iBasri et al.l 1 2010ll and 
iBasri et al ] (l201lh . However, because the Kepler pipeline was not designed to measure 
long-period stellar variability they used raw data from only a 90 day quarter, which 
introduced important limitations in their analysis. iBasri et al.l (2011) found that the 
amplitude of variability was in general larger for stars that were clearly periodic than 
for those that were not. As expected, the largest group of non-periodic variables were 
giant stars. To improve on the Kepler pipeline, which was designed for transit detection 
and not stellar variability, several groups developed their ow n correction software (e.g. 
lOarcia et aLlboiH i lMcOuillan et al.ll201 j i lOarcia et al.l[2014h . With these new analysis 
tools it was possible to detect longer-than-a-month rotation periods as well as to study 
the magnetic variability of the stars. 

Based on broadban d photometric variability, a new proxy for stellar magnetic activ - 
ity has been developed ( Garcia et al.l[201ol : IChaplin et alll2(lllal : lMathur et al.ll2014bl l. 
The time evolution of the variance of the signal appears to be a good indicator of the 
magnetic activity cycle. It has been successfully tested on the Sun using VIRGO/ SPM 
data recorded during the l ast 16.5 years (see for more details iGarcia et al.ll201.'i 1. Us¬ 
ing the same methodology. I^athur et al.l ( 2014al 'l studied a subset of 22 Kepler F stars 
with rotation periods of less than 12 days for which p-mode oscillations were reported. 
Because o f the relation between the rotation per iod of the star and the period of the 
cycle ('e.g. lBohm-Vitens3l2007l : l.louve etH ]S), they were expecting t o see in some 
stars cycle-like variations such as the one observed bv iMetcalfe et ahl ( 20101 1 in the 
F8V star i Horlogii. Th at star exhibited a ma gnetic cycle period of 1.6 years for a rota¬ 
tion period of 8.5 davs. l^athur et al.l (l2014al l found that only a small fraction of stars 
showed a regular cycle-like behavior (see FigO and there were no clear trends between 
activity and rotation in the F-star sample. A more systematic study still needs to be 
done to confirm or refute the existence of a general relationship between stellar cycle 
period and stellar rotation period for these hot stars. For solar analogs recent analy¬ 
sis of Kepler data have confirmed an rotation-age relationship ( do Nascimento et al.1 

I2OI4II . 


2 Global Models of Stellar Structure, Rotation, and Magnetism 

In order to understand, characterize and describe the physical processes at the origin of 
stellar convection, rotation and magnetism, with the Sun being an archetype of solar¬ 
like stars, we must develop a self-consistent, ideally nonlinear, description of their 
dynamics. Depending on the time scale of interest and given the current computing 
resources such model can be either 1-D, 2-D or 3-D. In the following section we briefly 
summarize the recent progress made in modeling the multi-scale dynamics of stars and 
their success at explaining the observed properties discussed in §1. 
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Fig. 5 Evolution of the normalized averaged variance of the Kepler F star KIC 12009504 for 
which p-mode oscillations have been measured (in black) and smoothed by a b oxcar function 
(160 d ays wide) in blue. A regular cycle of around 800 days is observed (see iMathur et al.l 
l2014al . for further details). 


2.1 1-D stellar evolution and core-envelope coupling 

The angular momentum evolution of a star - its rotational history from its birth until 
its death has profound consequences for its structure, evolution, and dynamics. Rota- 
tionally induced mixing for low ma ss stars can induce mild mixing that alters surface 
abundances (see IPinsonneau i3 ll997l for a review), while for massi ve stars, which rotate 
much more rapidly, the impact of mixing is far more dramatic (see lMaeder and Mevnetl 
|2012| . for a recent review). Stellar rotation is also intimately connected with the origin 
and generation of stellar magnetic fields, and the two together shape wind outflows 
from stars. Hence it is fundamental to be able to characterize it. The overall picture of 
angular momentum evolution was discussed in section 1.3. Here we focus on modeling. 

Structural evolution in stars will naturally generate internal angular velocity gra¬ 
dients, and in some evolutionary stages these shears can be large. Mass loss from 
the surface will also carry angular momentum away from the star, providing another 
powerful driver for radial differential rotation. In convective regions the interaction of 
rotation, magnetism, and convection will enforce an angular momentum distribution 
on a convective overturn timescale, much shorter than the natural nuclear evolution 
timescale. The response of radiatively stable regions is more complex, both because 
of our limited knowledge of the physics and because of limited empirical constraints 
on the relevant time scales. Internal redistribution by hydrodynamic and magnetohy¬ 
drodynamic physical process can both be important, and different classes of models 
testing different prescriptions have been employed. 
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Solar-like stars on the main sequence evolves on long time scales; the solar main 
sequence lifetime is of the order 10 Gyr. Computing such long temporal (secular) evolu¬ 
tion with multi-D codes is impractical. One must rely on 1-D stellar evolution models to 
describe th e structure, chemi cal and rotational evolution of stars. We refer to the recent 
review by ( Basil et al ][^ for a detailed discussion of 1-D stellar evolution codes, 
their success and difficulty at reproducing the solar internal structure, the influence of 
chemical composition, and how such helioseismically calibrated 1-D solar models have 
been extended to compute up-to-date 1-D stellar models. Here we focus on internal 
angular momentum transport and its consequences for the rotation history of stars. To 
this end, let’s apply the concept of angular momentum exchange and conservation to 
a star like our Sun. 

A useful fram ework that has been ext e nsivel y employed in the literature is the 
two zone model of lMacGregor and Bremi^ ( 199lh . In this model, the solar-like star is 
treated as a sphere divided into two zones: an inner stably stratified core and a turbulent 
convective envelope both independently rotating around an axis aligned with the Bz 
direction, where Bz is the unit vector along the z axis of a 3-D (x, y, z) cartesian system 
(see Fig[^left panel). Each possess their own moment of inertia I and angular velocity 
17; both structural evolution and magnetized winds can generate differential rotation. 
Their angular momentum is respectively Jcore — 7core l7core and Jenv — 7envl7env. 
Let’s now find the exact quantity of angular momentum exchange AJ between the two 
zones require to have them rotate uniformly at the rotation rate 17. The initial state 
is: = Jcore and = Jenv- Applying total angular momentum conservation e.g. 

= -^final + '^flnaP implies that the final state is: J^nal = -^core — AJ = /core 17 
and Jg„al -^env + AJ = 7envl7. A simple substitution yields: 


AJ = 


renv Jcore JcoreJf 


core<renv 


i + /e: 


( 6 ) 


If an external torque, such as the one coming from a stellar wind r^ind is applied to 
the star, the angular mome ntum evolution of the convect ive envelope and radiative 
interior can be expressed as ( MacGregor and Brennerlll99lh : 


JJcore AJ 

dt tc 

dJrmenv AJ Jenv 

dt tc tw 


(7) 

( 8 ) 


with tc the coupling time scale between the two zones due to the combined action of tur¬ 
bulence, waves, magnetic fields and viscous stresses and tw the wind braking timescale. 
These formula further assumes that the convective envelope transmits instantaneously 
the applied surface torque to the base of the convective envelope. 


Such models have been used to explore the relevant coupling time scale between the 
convective envelope and the radiative inter i or in solar-like stars over the course of their 


evolution ((MacGregor and Brennerl Il99ll: iKeooens et al. 1995: Krishnamurthi et al. 

1997 

lAllainI 

1998h and (for recent developments see Denissenkov et al.ll201fll: iBouvier 

2013 

iGallet and Bouvierll2013l:IOglethorDe and Garaudll2013l:IZhang and Penevl2014lL 


It is found that a time scale of tens to hundreds of Myr can explain the core-envelope 
coupling in young open cluster stars (see Figure O. 

One can identify key angular momentum transport processes which are expected 
to act inside a stellar convection zone or radiative interior: Reynolds and Maxwell 
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Fig. 6 Left: The two layers model of a stars, showing the exchange of angular momentum 
AJ between a core and an envelope rotating at different rates and the braking action of a 
stellar wind. Right: Evolutionnary track for t wo layer models showing how the coupling time 
scale can explain slow and fast rotators ifrom [Penissenkov et aO ll20inl ~l~l 


stresses, meridional circulation, viscous processes (usually negligible), large scale mag- 


netic toraue, gravity waves fsee Gilman et al.l 19891: IZahn 19921; Pinsonneault 

1997 


Mathis and ZahnI 2004. 200,4 Denissenkov et alJ 2010l: Mathid 2011 

: Bnin et al. 

2011 

Reiners and Mohantvl2012l: Mathis 2013bl;lGallet and Bouvieil2013l: 

van Saders and Pinsonneaultl 

2013|. and references thereinl. The amplitude and role of each process varies as a func- 


tion of time. For example, it is known that young stars are more active than old ones, 
so internal Maxwell stresses are expected to play a less significant role later in the star’s 
evolution on the main sequence and likewise for the stellar wind torque. Of course when 
the star undergoes a large structural change due to its secular evolution (as for instance 
when it goes up the red giant branch) the relative importance of different processes 
must be re-evaluated. 


iKawaletl 1 1988h proposed a genera l analytical fo rmulation of the applied torque that 
has been widely used. He adopted the jMes^ 196d}_^amework and adopted parameters 
tuned to reproduce Skumanich’s law ( Skumanich R97 for the time dependence of sur¬ 
face rotation. Many other studies have since improved upon Kawaler’s model by taking 
into account better observational information or theoretical improvements in order to 
assess how the magnetic torque applied by a stellar wind depends on stellar parameters 
and dynamics ( rotation, dynamo); see for e xample the Rossby-scaled saturation thresh¬ 
old adopted bv iKrishnamurthi et al.l (1993), which has been extensively employed sub¬ 
sequently. There ha s been much recent activity in this area, with new frameworks 
being proposed by (|Cranmer_aJid_Saaj 201l]; MahLgLuL 201j Reiners a nd Mohantvl 


201 Jivan Saders and Pinsonneaui^l201? BrowiJ 2014 : Matt~et al.l 201^ 


Self-consistent solar models including angular momentum transport from hydrody- _ 

nami c mechanisms were first generated in the 1980s ( Endal and Sofiall98ll : IPinsonneauIt et al.l 
Il989l 'l. Meridional circulation, hydrodynamic instabilities and gravity waves are the 
primary mixing and angular momentum tra nsport agents, so these models play a n 
important role in rotational mixing studies ( Mathij l2013al : ferun and Mathij l2014l l. 

Such models are sometimes co upled with ap proximations of magnetic effects, such as 
the Taylor-Spruit mechanism I Soruitl [2002l ~l . Thi s approach has been used to model 
mass i ve rotating stars ov er secular time scales ( Palacios et al.l [20031 : iMaraues et al.l 
I 2 OI 3 I : iMeibom et aI.ll2oii) and to generate grids of stellar models used for isochrones 
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and c hemical evolution studies I Ekstrom et al ] l2012l : lOeorev et al.ll201^ : iLagarde et alJ 
I 2 OI 2 I I. Some attempts have also been made to model in 2-D the secular evolution of 
internal angular mome ntum by considering the role of large scale magnetic field in the 
radiative interior fsee I Charbonneau and MacGregorl 1 19931 : iRuediger and KitchatinovI 
ll996l : [SDada et ^l201 fl and references therein). However such models do not solve for 
dynamo action as we will now discuss in the next section. 


2.2 2-D mean field stellar dynamo models 

On time scales of centuries to millenia it is possible to model in 2-D the magnetohy¬ 
drodynamics of stars using the mean field dynamo framework by parameterizing or 
neglecting some of the detailed description of convective motions. This allows one to 
explore the underlying physics and a wide range of stellar parameter space which we 
describe below. 


2.2.1 The Theoretical Framework 

Stellar magnetic cycles involve the generation and conversion of two components of the 
Sun’s magnetic field - the toroidal component and the poloidal component. Under the 
assumption of axisymmetry, the stellar magnetic field and plasma flow can be expressed 
as 


B — -f V X (^le^) 

V = r sm{d)n + Vp 


(9) 

( 10 ) 


In this axisymmetric decomposition, the first term on the R.H.S. of Eqn. 9 is the 
toroidal component and the second term the poloidal component of the magnetic field 
(of which A is the vector potential). The first term on the R.H.S. of Eqn. 10 signifies the 
^-component of the velocity (i.e., angular velocity 12), while Vp denotes the meridional 
{r-9) component of the velocity. 

Our understanding of solar and stellar cycles is base d primarily on 2-D kinematic 
0:12 dynamo models and full MHD s imulations (cf. S2 .2, |Brui^_et^ ] 1 I 2 OI 3 II (also pub¬ 
lished as chapter 5 of the ISSI book ( Thompson et allH and IWeisslll994l : ICharbonneaul 
1201011 ■ most of which have been constructed in the solar context. Here we briefly re¬ 
view this understanding. The toroidal component of the magnetic field is produced 
by stretching of an initial poloidal field by the differential rotation inside stars (the 
dynamo 12-effect) and is thought to be stored in a layer stable to convection near the 
base of the convection zone. Strong toroidal flux tubes rise up due to magnetic buoy¬ 
ancy and erupt as magnetic (star) spots. The toroidal flux must be recycled to poloidal 
flux for the dynamo cycle to perpetuate and this conversion mechanism (traditionally 
termed as the dynamo a-effect) could be achieved by various means. The two most 
popular an d widely stud ied mechanisms are the mean-field a-effect driven by helical 
turbulence I Parker Il955h and the decay and dispersal of tilted bipolar a ctive regions 
- i.e., bipolar sunspot or starspot pairs ( Babcock!Il96ll : lLeightonlll969h . hereby the 
BL mechanism. It is thought that strong toroidal flux tubes, believed to be necessary 
to produce sunspots, would quench the mean-field a-effect, if there is no segregation 
of the poloidal and toroidal generation process and thus the BL mechanism must be 
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the dominant mechanism for poloidal field generation in the Sun; however, this de¬ 
bate is far from over and we will return to this later. We display in Figure [7] a typical 
representation of a BL like dynamo model. 



Fig. 7 Left: a 2.5D kinematic dynamo simulation depicting a snapshot of the toroidal (right 
meridional plane) and poloidal field (left meridional place) configurations at a phase close to 
solar minimum. The simulation shows that the toroidal component of the field is generated 
and stored deep in the convection zone, from where it buoyantly erupts to form sunspots. 
The poloidal field is cre ated in the near-sur face layers in this particular dynamo model of the 
Babcock-Leighton type dNandv et aLllToilh . 3-D rendering of field lines (magenta) show the 
expected wound-up field structure in the solar interior. Right: Butterfly diagram for 2 models 
representative of pumping dominated dynamo, rotating at respectively 0.7f?Q and S.Of?©. We 
note the equatorward branch, the correct phase relationship between the toroidal and poloidal 
component of field (shown as vertical lines) and the sh orter rotation period for the case rotating 
at three time the solar rate llDo Cap and Brunll20llh . 


In the kinematic approach, the plasma velocity is prescribed and the magnetic field 
evolution equation is solved under the assumption of no direct magnetic feedback on 
the flows. The magnetic induction equation is given by 

^ = V X (v X B-? 7 V X B). (11) 


On substituting Eqn. 9 and Eqn. 10 in the induction equation and separating out the 
toroidal and poloidal magnetic field components we get the evolution equations for the 
poloidal and toroidal components, respectively: 


dB^ 

dt 


+ s 


Vp ■ V 


dA 1 , 


V(sA)] = rj (v^ 
+ (V • Vp)B0 = ?7 ( 


= r7 V" 


^)A-f5'„, (12) 

^) B^ + s{[V x(Ae^)] -VQ) 


1 d{sB^‘) dr) 
s dr dr 


( 13 ) 
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Here s = rsm{9) and rj is the turbulent diffusion coefficient. The source term for the 
poloidal field Sa could denote either the mean-field a-effect (in which case Sa = 
or the BL source due to the buoyant eruption and flux dispersal of tilted active regions. 
The flow in the meridional plane Vp represents meridional circulation (vp = Vm, where 
Vm denotes the meridional circulation), and in the presence of turbulent pumping 
(say v-y) it represents the combined effect of both (i.e., Vp = Vm + v-y). Appropriate 
prescriptions for the flows, turbulent diffusion and nature and profile of the poloidal 
field source along with the boundary conditions, completely describes the kinematic 
dynamo system defined in Eqn. 12 and Eqn. 13. 


In BL flux transport dynamo models, the cycle period Pcyc is found to depend 
strongly on the meridional circulation amplitude uq and its profile and less on the rota¬ 
tion r ate Qn or the amplitud e of the surface source term sg ( Dikpati and CharbonnearJ 
I 1999 I : I.Touve and BrunlboOTli : 


eye 


o0.05e0.07 -0.83 
OC fig >50 ^^g 


(14) 


As will be seen in the next section, the meridional circulation is found to decrease 
with the rotation rate as ug oc This is not intuitive as one could expect the 

meridional circulation velocity to increases with the rotation rate. A careful study of 
the 4> component of the vorticity equation shows that the longitudinal vorticity actually 
decreases with rotation rate as more and more kinetic energy is being transferred to 
longitudinal motions at the expense of meridional motions. The fact t hat in recent 3-D 
simulations ( Ballot et al.l I 2 OO 7 I : iBrown et al.l l2008l : iMatt et al.l I 2 OI 1 I : lAueustson et al] 
l2012l i. the meridional circulation is found to weaken as the models is rotated faster 
implies that standard advection dominated flux transport dynamo models yield the 
opposite dependency with rotation th an the one observe d, e.g. activity cycles are found 
to be longer for faster rotating stars ( Jouve et al jlimO). This fact alone requires us to 
revise our current dynamo paradigm for solar-like stars. One way is to short-circuit the 
advection path by, for in stance, adding more cells in latitude or increasing the radial 
diffusion as was done in | Jouve et al.l 2010|^ Hazraet_^ 2014|); another is to introduce 


magnetic turbulent pumping | Mason et al. 20081 : Guerr^o and de Gouveia Dal Pin3 


I 2 OO 8 I : [Po Gao and Brur]l201lh . These large-scale magnetic flux transport processes, 
e.g. meridional circulation, turbulent diffusion and turbulent (or topological) pump¬ 
ing of magnetic flux are extremely important for dynamo operation because they 
link the spatially segregated source layers of the dynamo, specifically coupling the 


2007; 


2012 ; 


influences the overall behavior of the dvnamo significantlv (I.Touve and Brun 

Yeates et al.l 2008 

; Guerrero and de Gouveia Dal Pinoll2008l: iKarak and Nandv 

Do Cao and Brun 

|2011^. Meridional circulation is also thought to critically in- 


fluence the nature of magnetic butterfly diagrams, including their structure, overlap be¬ 
tween cycles and cycle period s |cfliMbonneaiL_an(iJ3ikBatj|2000|; Nandv and Ghoudhuril 


I 2 OO 2 I ; Ijouve and Brunl|2007l : iNandv et al. 201ll: 


Hazra et al 


I 2 OI 4 I) and also impacts 


scaling relationships in stellar activity I Nandv 12004 ; Jouveet*^ 201(ih . However, the 
exact profile of the meridional flow remains unconstrained in star s. For the Sun, ev¬ 
idence for more than one meridional flow cell are being reported ( Haber et al ] l2002l : 
iMitra-Kraev and Thompsonll200'^ ; IZhao et al.ll2013h and there is serious debate about 
its relative role vis-a-vis other flux transport processes. 
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Fo r instance new stellar dynamo models including magnetic mimping lDo Cao and Brmi] 
(l 201 lh yield the following dependency for the cycle period Pcyc'- 

73 —0.40 —0.30 —0.15 /T r\ 

Pcyc OC Vg 7^q 7gQ (15) 

where jrO o-nd 7 eo represent the amplitude of the radial and latitudinal turbulent 
pumping terms. We see that the cycle period dependency on the meridional circulation 
is reduced and that if we want to recover the observational trends, in particular the 
shorter cycle period (again, we assume that iiq oc and assuming that 7rl7e 

remains constant, the pumping effect should roughly scale as Qg. Such a scaling may 
be too extreme and only systematic 3-D numerical simulations will tell us if this is the 
case or not. Nevertheless, pumping dominated stellar dynamos are a plausible solution 
to explain stellar observations. 


2.2.2 Stellar Activity and Magnetism: Confronting Mean Field Dynamo Models with 
Observations 


The general trends in stellar activity observations, in particular the rotation, cycle 
amplitude, and cycle period relationship is succinctly expressed when confronted with 
the stellar Rossby number, Ro = Pq /tc the ratio of the rotat i on period, Pg to the 
convective turn-over-time of the star jPurnev and Latouilll978l : iBaliunas et al.l [19961 : 
iMontesinos et al.ll200ll : Ijouve et al.ll201of )7 This number can be calculated using a com¬ 
bination of observations with some theoretical ideas about the nature of stellar convec¬ 
tion zones, assuming that the convective turn-over-time is Tc = LjV = /r], where L 

is a characteristic length scale, V a typical convective speed and 77 the magnetic dif- 
fusivity. The latter equality is valid if the magnetic Reynolds number Rm = VL/rj is 
taken to be of order unity, an assumption that depends on the properties of stellar con¬ 
vection zones and on the convection zone depth. It can be interpreted as characterizing 
turbulent mixing of magnetic fields in stellar interiors, while the rotation period plays a 
role in determining differential rotation. It is well known (both theoretically and obser- 
vationally) that the differential rotation in the con vective envelope of solar-type stars 


2005: Ballot et al. 

2007l:lBrown et al.ll2008l: iKiiker and Riidigei 

20081: Matt et al. 2011: 

Augustson et al. 2012h. However, the exact scaling m (i.e. AD 

oc Hg’") is still a matter 


of debate among the observers and theoreticians, being sensitive to the observational 
techniques used and to the modelling approach (see §2.2). Note that the o-effect used 
in mean field theory is related to helical turbulence and represents a measure or pa¬ 
rame terization of the electromotive mean force (emf) (e.g. lMoffattlll973 : IPouauet et al.l 
Il976l l. It is thus also linked to the rotation rate of the star and the amount of kinetic 
helicity present in its convective envelope. However in Babcock-Leighton flux transport 
dynamo models the classical a-effect is replaced by a surface term linked to the tilt of 
the active regions with respect to the east-west direction (the so-called Joy’s law, e.p 


I Hale et lBabcocl3ll96ll: iKosovichev and Stenflolf2008l: IPasi-Espuig et ^l 201 

This tilt is thought to be due to the action of t he Coriolis force during the ri se of the 
toroidal structures that emerge as active regions JP’Silva and Choudhurilll993lL Recent 
3-D simulations in sph erical shells Jjouve^n^^ruiJ^ 0^ Fanll200^ Weber et al.ll201ll : 
IPinto and Bnmll2ni3l : l.louve et al.ll 2 ni^ Weber^t al.l 2013l ~) seem to indicate that this 
is not the only effect responsible for the observed tilt and that the twist and arching 
of the toroidal structures as well as the continuous action of the surface convection 
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during the emergence influence the resulting tilt. The Rossby number, overall, char¬ 
acterizes stellar convection zones - the zone where the dynamo mechanism generates 
stellar magnetic fields. Not surprisingly, the Rossby number is simply related through 
scaling arguments to the dynamo number, D ~ IjRo^. Thus the Rossby nu mber can be 
utilized to relate stellar observations to the dynamo mechanism (see e.g.. iNoves et al.) 
ll984bl[T984al : lTobiasiri998l : iMontesinos et ahlboOll : lNandvir2004h . 

On can also implement a non linea r feed back via the l arge s cale Lorentz force 
(also called the Malkus-Proctor effect; iMalkus and Proctor! JlfiTSh l on the dynamo 
field generation. This has the advantage to y ield a satura t ion (quenching) process 
based on the actual dynamics of the models ( Tobiad 1 19971 : IMoss and Brool3 I 2 OO 0 I : 
lBushbvll2006l . and references therein). These models demonstrate that for low magnetic 
Prandtl number (e.g. Pm = v/r]), such that the viscous time scale is much longer 
than the magnetic diffusion time scale, the dynamo is irregular, with grand minimum¬ 
like behavior and switching of parity between symm etric and antisymmetric dynamo 
families as observed in the Sun dPeRosa et al.ll2012l l. 

From the observational perspective, analysis shows that stars with low Rossby 
numbers typically exhibit strong activity amplitudes and have irregular cycles, while 
those stars with relatively higher Rossby numbers tend to have lower activity levels, are 
more likely to host magnetic cycles and are more likely to be found in grand minima like 
phases. Since lower Rossby number indicates a higher, underlying dynamo number, and 
vice-versa, this result is consistent with dynamo theory. Stars which host very efficient 
dynamos (having high P) are expected to have high activity levels; the non-linear 
nature of the dynamo mechanism, in conjunction with a high dynamo number also 
generates irregular cycles because the magnetic feedback is likely stronger. On the 
other hand, stars with low or moderate dynamo numbers are relatively less active, and 
are more likely to host regular cyclic behavior with occasional Maunder-like states - 
such as the Sun. The consistency of stellar observations with the theoretically expected 
behaviour of a MHD dynamo mechanism leaves little doubt that internal dynamos are 
the ultimate source of stellar activity. 

A closer look at the activity period versus rotation rate relationship yields more 
stringent constraints on the nature of stellar dynamos. In particular, it is seen that the 
periods of stellar magnetic cycles decrease with increasing rotation rate (i.e., decreas¬ 
ing rotation period Pfi). As discussed in §2.1, there is an inter-dependance between 
dynamo action, magnetic activity levels, rotation and stellar evolution via the torque 
that magnetized thermally driven stellar winds apply on s olar-like stars. The mass de- 
pendance of this torque is subject to intense studies (see ISouviei] l2013l . for a recent 
discussion). It is understood that solar dynamo models that rely on the BL mechanism 
for poloidal field generation are critically dependent on flux transport processes such 
as meridional circulation which couples the source regions for the toroidal and poloidal 
components of the magnetic field. Such BL dynamo models have been successful in 
reproducing various features of the solar cycle. Based on such a BL dynamo model, 
iNandvl 1 2004l j has shown that the stellar activity amplitude and period versus rotation 
rate relationship can only be reproduced if both the rotational shear and meridional 
flow speed scales positively with rotation rate. While full MHD numerical simulations 
show that the former holds, t hey also indicate that meridional circulation spee d does 
not increase with rotate rate ( Augustson et al ] l2012ll . IDo Cap and Bruiil 1 I 2 OI 1 II point 
out in a recent study that the inclusion of turbulent pumping, with the speed of the 
former increasing with rotation rate, is necessary to recover t he stellar scaling laws in 
the light of these new developments in full MHD simulations. IPube and Charbonneaul 
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have also computed 2-D kinematic dynamo using ingredients (such as the alpha 
effect) deduced from 3-D simulations. They show that essential ingredients of the com¬ 
plex 3-D dynamics (such as for instance the cycle period or the sense of propagation 
of the dynamo wave) can be recovered via such techniques. 3-D k inematic dynamo 
models or convective B-L models have also recently been developed ( Chan et al ] l2007l : 
iMiesch and Brownll2012| : iMiesch and DikDatill20l3l . 

Hence confronting stellar activity observations with diverse simulation techniques 
is very useful and allows the community to test new ideas and concepts. 

We now discuss in more details nonlinear, 3-D MHD simulations of solar convection 
zone and magnetism. 


2.3 3-D Global Models of Stellar Convection and Dynamo Action 


With the advent of massively parallel computers it is now possible to attack the diffi¬ 
cult problem of stellar dynamics and magnetism with direct 3-D nonlinear numerical 
simulations. As is frequently the case with multidimensional MHD numerical simula¬ 
tions, two approaches are possible: local high re solution simulat i ons that describe the 
smal l scale turbulent motions and local dynamo ( Cat,tan^ll999l : |Vogler and Schiisslerl 
I 2 OO 7 I I and global ones that take into account t he correct geometry and t he existence 
of larg e scale flows in rotating con vective zones iGlatzmaieil ( 1985alll985bl ') ; iBrun et al.l 
( 2004 1: IMiesch and Toomrel 1 I 2 OO 9 I I. Given the broadness of the topic, here we present 
the progress made in modeling stellar convection and dynamos with the global ap¬ 
proach. These simulations allow us to describe turbulent, rotating and magnetized 
convection either in a full sphere, a shell or a wedge-like geometry (cf. the recent 


book by iGharbonneaul l2013al l . Gontrary to the kinematic mean field stellar dynamo 
approach discussed in the previous sections, 3-D simulations use the full set of MHD 
equations (we refer to these recent references for simpler 2-D mo dels of stellar differ¬ 
entia l rotation assuming a parametrization of turbulent convection Kiiker and Riidigerl 
I 2 OO 8 I : iHotta and Yokovamal I 2 OI 1 I : iKitchatinov and Olemsk^ I 2 OI 2 I 1 . These equations 
couple classical fluid mechanic equations to the induction equatio n and as such do 
not prescribe the velocity field but compute it self consistently (c.f. iBrun et~S^ I l2004 
Im ieschl [ 200 ^ . This is the main advantage of the global full MHD approach: the in¬ 
fluence of rotation and magnetic field on turbulence and convection is computed not 
parametrized. 

Over the last three decades several groups have developed stellar convection and 
dynamo models in spherical geometr y following the pioneerin g work of P. Gilman and 


try t( 


G. G l atzmaier in the early-mid 8 0’s llGilman an d Milleiiri98ll : iGlatzmaier and GilmanI 


1 19821 : iGilmanll 19831 : lGIatzmaieijll985a 1985^1 1" Prominent examples of modern global or 
wedge-like parallel codes that haye in recent years published original studies on stellar 
conyection and dynamo include: 


the Anelastic Spherical Harmonics (ASH'! code: IMiesch et al. (1200 

"): Brun and Toomrel 

(2002 

hiBallot et al.l (l200itl: Miesch et al.l(20081'): Browning! (I 2 OO 8 I') 

Bessolaz and Brun! 

(2011 

1: Augustson et al. 

(20121 on surface convection and zonal flows and iBrurJ 

(2004 

1: Brun et al. 

(2004 

1: Browning et al. ( 

200fill: Brown et al. ( 

2008 ). 2 OIOI. 


Brun et al.l (201T): 

Miesch and BrowrJ (I2012 

hiNelson et al.l (120131 

:lAugustson et al.l 

( 2013f): iNelson et al.l ( 2014h on solar-like star dynamo action, cycles and magnetic 


wreaths, 
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— the Eulag MHD code: lOhizaru et alJ (201 fll ~): iRacine et alJ (l201ll ~): ICharbonneaij 

1 2013bh on solar an d stellar dynamo and cy c les, _ 

— the PENCIL code: IWarnecke et alJ 1 2013li : iKaovla et alJ ( 2013l l: ICuerrero et alJ 
( 2013ll on stellar convection, dynamos and cycl es, 


— the Magics code: 


Christensen and Aubert| 1 200fill : [christensen et aO ( 2009h : IChristensenl 


1 2O10h : lGastine e taLTif^ oT^ 20141 on convection and dynamos 


— K. Chan’s code Cai et alJ (1201 111 on core convection 


Other codes have been mostly used to model the geodynamo with some applications 
to the dynamo in M dwarfs. 


— Leeds code 


byE 


— Glatzmaier’s code: adapte d to model Earth’s dynam o lClatzmaier and RobertsI (19^) 
or giant pl anets dynam ics Glatzmaier et all 1 20091 ) . 

— EBE code: [chajj_et_al. (|2007|), 


H^eidraeh and .lonesi ( IDfl.lh : [Sreenivasan and .TonesI 1 20061 ') ; dynamo 
id ll201lll 


— PARODY code: ICoudard and Dori^ I 2008h : iMorin et al. | 201lh: jSchrinner et al.l 
( 2012I . I 2 OI 4 I I on geo and stellar dynamos and Aubert_et_^ (12007 ) on the role of 
thermal wind in setting me an large scale flow s and Aubert_eLalJ 1200^ o n dyn amo 

— Busse and Simitev’s code: ICrote and Bussel ( 2001 1: ISimitev and Bussd J200illl on 

dynamo and bistability __ 

— Anrnou and Heimpel’s code: lHeimpel et al.l 1 200^ : lAurnou et al.l (20o3) on dynamo 

in convective envelopes __ 

— Severa l research groups in Ja pan: lTakehiro and Havashil ( lOOflI l: Yoshid^^n^j^agewmal 


Kirnur^^£_ay_([20n ) on conv e ction or zonal flows or iKaeevama and Satd 


1997 1: iNishikawaandl^sand ( 2008l l : iMasada et al.l ( 2013l l on dynamo, solar con¬ 
vection and magnetism, the most recent studies making nse of a novel “ying-yang” 


grid to map a spherical shell ( Kaeevama and Satdl20041 . 


Many of the codes nsed in the geophysics community were originally developed to model 
the geodynamo. They therefore rely on the Boussinesq approximation, which is not 
appropriate for highly stratified and compressible convective envelope as found in solar¬ 
like stars. Modelers have been converting these codes t o be able to d e al with stellar-like 
stratified flow by using the anelastic approximation ( GoughIIl969l : Iciatzmaieil Il984l : 
IVasUet^. 201T1 and validating their codes with a recent international benchmark 


( Jones et al. 2 OI 1 II . 


Given the intrinsic complexity of computing nonlinear convection and dynamo sim¬ 
ulations and the large diversity of stars to model, we first b riefly describe pu rely hy- 
drodynamical models (see also chapter 3 in the ISSI book ( Thompson et akh i. Such 
models reveal the fundamental ingredients at the origin of large scale mean flow, such 
as differential rotation and meridional circulation in rotating global convection models. 
We then follow up by describing the origin of magnetic fields (either cyclic or irregular) 
through nonlinear dynamo action and the feedback of such fields on the flows - making 
the link with the 2-D models discussed earlier when necessary. 


2.3.1 Convection and Large Scale Flows 

Convection is present in stars of all masses on the main sequence. It takes either 
the form of a convective core in massive ones (greater than about 1.3 Mq) or of a 
convective envelope of varying thickness for lows mass stars (from early F down to 
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about M3), with stars less massive than about 0.3 Mq being fully convective. Here we 
will limit our discussion to solar-like stars possessing a deep convective envelope. For a 
discussion on two other categories (e.g. high mass and very low mass stars) interested 


readers may refer to (e.g. Browning et al. 

to 

o 

o 

Brun et al.l 

2005: Dobler et al. 200fil: 

Brown et al. 20081: Christensen et al. 200f 

: Featherstone et al.l 20091: Augustson et al.l 

201ll: Cai et al. 201ll:lMorin et al. 2011: ISchrinner et al. 20121. and references therein). 




Fig. 8 Stellar luminosity, effective temperature and mass contained in the convective envelope 
vs spectral type in solar like stars with the solar metallicity (computed with the CESAM code 
lMore]lll997l'l 



Fig. 9 3-D rendering of turbulen t surface convection in solar-like stars with varying convection 
depths llBessolaz and BrunI 1201111 . The radial velocity is shown with cold downflows in blue 
tones. 


Convection plays a key role in transporting the star’s heat by having strong corre¬ 
lations between vertical motions and temperature fluctuations. It is a fundamentally 
turbulent process in stars as the Rayleigh number {Ra, e.g. the ratio between buoy¬ 
ancy driving and dissipation process) is huge, of order 10^® in the Sun. This leads to 
a large range of scales from possibly giant cells on the order of the thickness of the 
convective layer (e..g 200 Mm in the Sun) down to granulation size (1 Mm), with in¬ 
termediate scales s uch as supergranulation (20 -50 Mm) well delineated by the surface 
magnetic network i Rieutord and Rinconll^ioh . As the luminosity (energy output) of 
the star changes so does the amplitude of the convective flow. A simple scaling argu¬ 
ment based on mixing length treatment (MLT) gives Vr oc [L^:/{pczR^))^^^, with L* 
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and R* respectively t he star’s luminosity and radius and pcz the averaged density of the 
convective envelope ( Kipnenhahn et ^l2013l l. Since the luminosity varies by at least 
2 orders of magnitude across the spectral type we are considering (i.e early F to early 
M) {L^, oc Mf), the stellar radius only by a few (i?* oc M®'®), it is expected that the 
vigor of the convective flows do likewise modulo the effect of density. As can been seen 
in Figure 8, 1-D stellar structure models of solar-like stars indicate that the averaged 
mean density changes significantly in the convective envelope of these stars, in such 
a wa y that more mass ive stars have on average a lower mean density than low mass 
ones ( Matt et al ]H). This is easily understood as the base of the envelope moves up 
towards the surface with the thinning of the convective envelope. This yields a large 
amplitude difference between F and K stars for their characteristic convection velocity. 
This difference in velocity amplitude is recovered in global numerical simulations of 
stellar convection. We illustrate an example of such 3-D turbulent convective simula¬ 
tions in Figure [9l where we clearly see the complex network of downflows surrounding 
broad upflows. 

The turbulent motions within the convective envelopes of rotating stars experience 
the effect of the Coriolis force to a degree that depends on their fluid Rossby number 
Ro = here define d as the ratio between turbulent and the so-called planetary 

vorticity ( Pedlosk'^l 1982h . We note that several deflnitions of the Rossby number exist 
in the community (stellar Rossby based on the ratio between convective turnover time 
and stellar rotation period; convective Rossby-^ Ra/TaPv, local Rossby number), but 
it can be shown that they depend almost linearly on one anot her (see for example 
iKim and Demarauel Il996l : iLandin et al.l I 2 OIC 1 I : iMatt et al.l I 2 OI 4 I I . We here choose to 
use the fluid Rossby Ro because a clear transition of behavior occurs at Ro ~ 1 
with the differential rotation becoming prograde or retrograde. Indeed if the stars 
rotate fast enough such that the large scale motions (larger than the Rossby radius of 
deformation) feel the continuous action of the Coriolis force as they rise and sink then 
the redistribution of angular momentum is such that a non uniform rotation profile is 
established in the convective envelope. 

In figure [TO] we display 4 realizations of the longitudinally av eraged differentia l 
rotation and meridional circulation computed with the ASH code ( Matt et al.ll2014l . 
in preparation). These profiles have been obtained by changing either the stellar mass 
(i.e aspect ratio and luminosity) or the reference frame rotation rate. We note several 
trends: As the rotation rate is increased, the differential rotation becomes more and 
more aligned with the rotation axis (along Cz) and likewise the meridional circulation 
posses more and more cells both in latitude and radius. The 1.1 solar mass case rotating 
at one time solar rate exhibits an anti-solar differential rotation, with slow equator 
and fast pole. It is indeed the only case with a Rossby number greater than 1 as 
seen i n many recent studi es jBessolaz and Brurj r201ll : iMatt et al.ll201l] : ICuerrero et al.l 
l2013l : [Gastine et al.ll2ni'i : lKaDvla et al.l 2014h . We thus see that the transition between 
prograde and retrograde differential rotation will not occur at the same rotation rate 
(v sin i) for each stellar spectral type. We also note the strong shear l ayer at the base 
of the convective domain (e.g. the tachocline ISniegel and Zahnlll99^ that naturally 
develops in each modeled stars with a profile closer to the solar case for all the cases 
with a Rossby number Ro ~ 0.6. 

These profiles can be understood by considering the transport of angular momen¬ 
tum i n a spherical convective shell. When perfor ming this analysis as in ( Brun and Toomra 
I 2 OO 2 I : iBrown et al.ll20oi : lAueustson et al.l|2012l ') , we find that the Reynolds stresses are 
at the origin of the prograde equator for models with a low Rossby number. Through 
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Fig. 10 Differential rotation and meridional circulation realized in 3-D models of solar-like 
convective envelop e coupled nonline arly (i.e both zones back react on one another) to a stably 
stratified interior (iMatt et al.l 12014 . in preparation)). 3 masses (0.5, 0.9 and 1.1) Mq and 3 
rotation rates (solar, 3 and 5 time solar) are being shown. Starting from the left in each row: 
color contours showing Q{r, 9) and the meridional streamfunction for a 0.5 Mq star rotating at 
5 Hq, 0.9 Mq star rotating at 3 and 1 Hq and 1.1 Mg star rotating at 17q. Prograde rotation 
is shown in red/white colors and retrograde rotation in dark/blue. Clockwise meridional flows 
in red. 


gyroscopic pumping the meridional circulation transport angular momentum to coun- 


the contribution from viscous effect). We refer to (ZahrJ 

19921: M 

athis and Zahn 2004: 

McIntvrel2007l:lGaraud and BodenheimeJ2010l:lBrun et al.ll201ll: 

Miesch and HindmanI 

2011 ) for a more thorough discussion of the role of the meridional circulation in reach- 


ing angular momentum balance in rotating stars. The orientation of the iso-contours 
of 17 from being cylindrical to being more conical (solar-like) at mid latitude is due 
to the effect of an efficient thermal wind associa ted with latitudinal entropy variations 
(with relatively hotter poles and cooler equator) [Miescl^_et_alJj200^)Trhese variations 
are enhanced further by the presence of a tachocline llBruiTetaLl^llI ). 


One important aspect of performing 3-D numerical simulations is that one can ob¬ 
tain general trends on the variations of the main physical processes acting in convective 
envelope with various parameters such as rotation rate, luminosity, aspect ratio. We 
focus on scaling laws for the differential rotation and meridional circulation by synthe¬ 
sizing the recently published 3-D stellar convection models computed with the ASH 
code. The models rang e from 0.5 up to 1.3 times the solar mass and from 1 to 20 ti mes 
the solar rotation rate ( Brown et al.l[200^ : iMatt et al.l[201ll : lAueustson et al.ll2oi3) . In 
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Fig. 11 Differential rotation Afl as a fun ction of Rossby nu mber for stellar masses ranging 
from 0.5 up to 1.3 1 cases from dBrown et al.|[2b08n . 1.2 and 1.3 Mq cases coming 

from llAuEfustson et al]l2012f) V Kinetic energy of the axisymmetric meridion al circulation as a 
function of the Rossby number for models ranging from 0.5 to 1.1 Mq, from JMatt et al.ll20l3 . 
in preparation). 


Figure[TT]we show how the latitudinal differential rotation contrast varies as a function 
of the Rossby number. There is a clear trend for AO to rise with increasing rotation 
rate. This is in qualitative agreement with direct and indirect observational data that 
imply larger absolute shear rates in rapid rotators to produce stronger magnetic fields 
and enhanced angular mome ntum loss relative t o slow rotators. The actual slope seems 
a bit larger t han observed by Collier CameronI 1 I 2 OO 7 II but not so when compared with 
the work by IPonahue et al.l ( 19961) . SaaiTl^llh found that the presence of binaries 
in the published samples, which can be distinct from single stars because of tidal ef¬ 
fects have led to an underestimate of the dependency of AQ on rotation rate. Clearly 
more observational data are required to constrain better the exponent n in the equa¬ 
tion AS7 oc 17". We also clearly see that more massive stars possess larger ASl for a 
gi ven Rossby number. This is in reasonable agreement with the observations published 
in iBarnes et al.l (20^). Although, given the paucity of data and the rather large error 
bars, the exact exponent is uncertain. Data from the K epler mission and in the (future) 
PLAT O mission will be a major asset in this context I Garcia et al ] |2014l : iRauer et ahl 

I 2 OI 3 I) . 


Another interesting trait is that the meridional circulation is found to decrease 
in strength as the stars are made to rotate faster. Rapid rotation causes more kinetic 
energy to go into toroidal motions, which can be anticipated by looki ng at the equation 
for gyroscopic pumping ( Landin et al ] I 2 OI 0 I : lAugustson et ahl I 2 OI 2 I ) . This has a very 
important consequence: Flux-transport dynamo models as discussed in the previous 
section can’t be easily extended to other stellar spectral type to explain the dependence 
of their cycle period with their rotation period. Indeed flux-transport models of the 
advective type depends directly on the meridional circulation amplitude to set the cycle 
period. Observations indicate that the cycle period should decrease with increasing 
rotation rate. However, given that the meridional flow is fo und to weaken, th is can 
only result in an increase of the activity cycle as was shown bv l.Touve et ahl ( 20101) . One 
needs either to tak e into account multi-cellular flows as obser ved also in 3-D simulation s 
( Jouve et al.ll201ol) or take into account turbulent pumping I Do Cap and BruDll201ll ). 
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2.3.2 Nonlinear Dynamo Effect, Magnetic Activity, and Cycles 


The observed magnetic activity shows the existence and period of cycles seem to depend 
on stellar spectral type and rotation rate. An obvious trend is that old (older than the 
age of the Hyades cluster ~ 800 Myr), slowly rotating solar-like stars are less active 
than their younger counterparts. We have also discussed in the previous section that 
convection and large scale flows change with spectral type. It is certainly not surprising 
that the associated dynamo action and magnetic activity do so as well. 

Several groups have recently published extensive parameter studies of dynamo ac¬ 
tion in rotating convective spherical shells of various thickness and degree of strat¬ 
ifications. The key results can be summarized as follow: 3-D numerical simulations 
of dynamo action in solar-like stars have revealed a large range of behavior , from 


steady dynamo, to i rregular and cyclic o nes |lBrurL_et_aL [2QQ4 


Racine et al. 

2011 

Kanvla et al. 

2013 


iGastine et al.l 


20 ^: _ „ . 

2015 ; [Neisoiret*^ 20131 : ISchriimeTmiL 2014 , and references therein). 


Schrinner et al.l 12012 


Brown et al. 


Augustson et al 


2O10l 

et al.l 


2011 


2013 


For low mass (M-type) stars, it is seen that a bistab i lity si milar to that found for 
the geodynamo is possible. Observations bv iMorin et al.l ( 2010li seems to indicate that 
a strong dipolar state and a multipolar state may coexist for stars having similar stellar 
parameters. Dynamo models show that around a local Ross by number of 0.1 a weak 
(mult ip olar) and strong (di polar) dynamo branch may coexist Ichristensen and AubertI 
(20^); iMorin et al.l ll201ll 'l. Such bistable states are found in a given range of pa¬ 
rameters (low stratification, high Pm) and may not actually exist in real stars. It 
may be the case that these stars are observed in a different phase (low vs high) of 
their activity cycle. Other similarity between the geodynamo and low mass stellar dy- 
namos also exist such as the scaling of the magnetic energy with the available heat flux 
( Christensen et alJ 12001^ 1. Such scaling tends to break down for higher mass stars for 
which the tachocline may play a significant role. 

Others 3-D simulations of F, G and K stars show that for high rotation rates, large 
scale magnetic wreaths (see Figure left panel) are obtained without requiring the 


presence of a tachocline iBrown et al. 


( 201Cll . [^llh . We show on Figure [12] four butterfly 


diagrams (time-latitude plots of the azimuthally averaged toroidal magnetic field near 
the base of the CZ) realized in such simu lations for of a solar-like star rotating at three 
times the solar rate [Nelson et all 1 2013li . We remark that as the model is made more 
turbulent, the steady magnetic wreaths become more time dependent and can lead to 
cyclic activity (bottom right panel). 

In such stars, along with the degree of turbulence, rotation plays an important 
role in determining the global properties of their magnetism. This is due to a shift 
in the balance of forces driving the flow between the advection, Coriolis and Lorentz 
terms. As the rotation rate increases the Lorentz force tends to balance the Coriolis 
force yielding larger magnetic energy in superequipartion with the kinetic energy of 
the flow (a direct consequence of a magnetostrophic state; c.f. strong scaling below) 
as in the Earth’s iron core. The Elsasser number A = /AtvPczVDq, with pcz mean 
density in the convective envelope, p magnetic diffusivity, i?o stellar rotation rate, B 
a characteristic magnetic field of the CZ, is useful to discuss this balance of terms in 
the Navier-Stokes (N.V.) equation. Depending on the amplitude of this number and 
on the balance assumed in the N avier-Stokes equation, various scaling of the magnetic 
field amplitude can be expected ( Fauve and Petrelisll2007l : IChristensenll201^ : 
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Fig. 12 Magnetic wreaths yielding in turn steady (case D3; top left), irregular (case D3a; top 
right) and quasi cy clic (cases D3b Sz S3; bottom left and right) magnetic butterfly diagrams 
iNelson et al.l ll2Q13l) . Red tones correspond to positive toroidal polarity. 



Fig. 13 Magnetic wreaths shown by 3-D rendering of magnetic fleld lines (iBrown et al.ll201Qr) . 
Temporal sequence of meridional cuts showing the mean poloidal field realized in case D5 of 
IBrown et al.l ||2011J) which was extrpolated using potentiel field approximation to about 1.5 
stellar radius. Note the global polarity field reversal of the poloidal field. 


— To an order of magnitude, an estimate of the equilibrium magnetic field (assuming 
ideal gas law) can easily be obtained: Beq ~ -^SnPgas ~ \Jpcz-, since T^ff varies 


by a factor 2 to 3 b etween earl; 
than a factor of 100 (Matt et al, 


F an d late K stars, whereas pcz varies by more 


y r an a . 


If one assumes for the sake of simplicity that the magnetic Reynolds number Rm = 
vd/rj ~1 such that a characteristic velocity is given by u ~ ?7/d0, the balance of terms 
in Navier-Stokes Equation is then given by: 

— Laminar (weak) scaling: Lorentz ~ viscous diffusion 

^ Bleak ~ PczW/d ~ pezV-nlS 

— Turbulent (equipartition) scaling: Lorentz ~ advection 

^ ^turb ^ PczV ~ pczp /d ^ l^turbl^rn 

— Magnetostrophic (strong) scaling (e.g. Elsasser nb yl ~ 1): Lorentz ~ Coriolis 

^ ^strong ^ pcz^Op 

^ note however that such choice of characteristic velocity scale is subject to caution in stellar 
turbulent convective envelopes that usually possess very large Rm 
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with V, d are characteristic velocity and length scales, Pm = ^Irt the magnetic Prandtl 
number. Of course there is an upper limit to the magnitude of the magnetic energy 
ultimately set by the amount of energy (likel y the star’s outward energy flux) that 
can b e made available to the dynamo process ( Schubert et al.l[l996l : IChristensen et ahl 
I 2 OO 9 I I. We recall here th at dynamo ac tion does not exist for all fluid motions due to its 
intrinsic 3-D character I Moffattlll978h . 

So a possible scenario for stellar dynamos based on heuristic scaling arguments 
could be the following: Stars rotating at moderate rate (such that their Elsasser number 
is small), have a level of magnetic energy (or averaged global field strength) that are 
less than or of the order of the equipartition field given by either the weak or turbulent 
scalings described above. As stars rotate faster and get closer to a magnetostrophic 
state with an Elssasser number of order 1 or larger, the formation of large and intense 
magnetic wreaths starts. The magnetic field becomes more and more dominated by its 
toroidal component and the magnetic energy becomes larger and larger going above 
the equipartition value set by the turbulent scaling as it now follows the strong scaling 
which depends linearly on the rotation rate. The consequence may be the following. As 
the magnetic energy (or field amplitude) becomes large, the associated Lorentz force 
starts back reacting strongly on the mean flow. The first consequence is what can be 
called an “f2-quenching”, e.g. the differential rotation reduces in strength and an a lmost 
solid body rotation sta te in the convective zone (envelope or core) is established I BrurJ 
l2004l : [Brun et allboosl i. Using mean-field classification, one may then hypothesize that 
the stellar dynamo transitions from an a — 17 or — 17 to being an o? dynamo, 
i.e. helical turbulence is solely responsible for field generation and maintenance, the 
large scale shear now plays a more marginal role. At that stage what remains of the 
magnetic wreaths is still unclear, more work must be done. As the rotation is made 
even faster, quenching of the a effect, due to the large scale magnetic field being more 
and more intense, occurs naturally as now the Lorentz force back-reacts directly on 
the helical convective motions and not solely on the large-scale mean flow. Following 
this discussion, one question naturally crops up which requires detailed investigations: 
How are these diverse set of plausible dynamo scaling relationships (involving magnetic 
activity amplitudes) related to stellar activity trends and can the latter impose more 
stringent constraints on solar and stellar dynamo processes? 

As we have seen in section §1.6, the activity of solar-like stars can be assessed by 
observing their X-ray luminosity Lx- Using such proxy it has been shown that magnetic 
activity increases with rotation rate and that above a given rotation rate (that depends 
on the st ellar spectral type but is around Ros ~ 0.1), stellar magnetic activity first 
saturates ( Pizzolato et al.ll2003l : IWright et al.ll201l[ l and then possibly “over-saturates” 
for extremely fast rotation rate s uch as the ones found in young solar-like T-Tauri stars 
( Feigelson and Montmerl^ IiOOOI L One can thus simply seek to characterize how the 
saturation of the dynamo generated field can be linked to the saturation of solar-like 
star X-ray emissivity. Such a link is not as straightforward to deduce as one could 
anticipate. Indeed one must also assess how the filling factor / of the magnetic field 
on the star’s surface evolves with stellar parameters not just the field strength. As we 
have noted the field amplitude does not easily saturate since the Lorentz force can 
achieve a dynamical balance with various terms whose importance and amplitude vary 
with stellar parameters. One may conjecture that the first saturation of the X-ray 
luminosity is due to f7-quenching and to the physical limitation of the spot coverage 
(e.g., filling factor /) on the stellar surface and that the second “over-saturati on” 
may be due to “a-quenching”, so of the actual field strength (see ICondoinl I 2 OI 2 I for 
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an alternative explanation). Currently no global-scale 3-D nonlinear MHD dynamos 
actually include spot formation and emergence. Global flux emergence simulations 
have been computed to assess the role of the Coriolis force, Hoop stresses, convective 
motions, large scale mean flows and background magnetic fields on the evolution and 
morph ology of simplified active regions by varying the properties of idealized global flux 


ropes JPanlboOtil: 


Jouve and Bru: 


re gions 

320091: 


Weber et al 


l e properties or idealized global tlux 
] l201ll : llsik et ahlboill : Ijouve et al.l 


I 2 OI 3 I : IPinto and BrurJl2013l : IWeber et al.l 2013h but until very recently, none had self- 
consistently dealt with the formation of sunspots in a cyclic convective dynamo. So in 
order to be able to physically characterize the saturation of B and / independently, 
spot-dynamos - wherein dynamos self consistently produce rising omega-loops - must 
be developed and the parameter space explored systematically. 

There has been recent progress in this direction using improved numerical treat¬ 
ment in 3-D simulations of the diffusivities (via dynamical Smagorinsky, slope lim¬ 
ited diffusion or implicit large eddy simulation (ILES)). Authors of these studies have 
been able to lower the diffusivity level low enough to yield more time dependent dy¬ 
namo solutions possessin g both cyclic behav i or, regular butterfly diagram and torsional 


namo solutions possessin g botn cyclic benav i or, regular butterny diagram and torsional 
oscillation-like behavior jKajada/ etjah 20n ^ 201^ 20^ DubeIn^_Qiarbonneau|20^; 


ISmolarkiewicz and CharbonneaJ 2013 ; Beaudoin et al.l 20131 : Augustson et_ah J2013|l. 

Some have even obtained for the first time buoyant magnetic wreaths ( Nelson et al.l 
I 2 OI 1 I . [ 2 OI 3 I I. In these recent simulations, magnetic wreath-like structures become tur¬ 
bulent and intermittent enough, that many intense bundles of fields reaching 50 kC 
start to form f7-loop like structures with statistical pro perties for rising loops that 
qualitatively agree with observations ( Nelson et al.ll2014l ~). We believe that such simu¬ 
lations are the progenitors of future more realistic spot-dynamos that will allow us to 
characterize the link between dynamo action, flux emergence and large scale magnetic 
topology as a function of stellar parameters. 

Finally also note that the torque applied by the magnetized wind (see discussion in 
§1.5 and 2.1) depends mostly on the low order fields and that the field to pology is ex¬ 
pected to change as a fun ction of stellar cycle and overall rotation rate ('see lPinto et al.l 
l2nill : iReville et aklboidh . Hence assessing the large-scale magnetic topology of stellar 
magnetic fields and how it changes as a function of stellar parameters is very important 
as it has a direct bearing on stellar evolution. 


3 Future Directions 


The characterization of the many facets of stellar dynamics (convection, rotation, mag¬ 
netism) has made tremendous progress over the last 20 years. We have moved away 
from a quasi-static view of stars to a time dependent, turbulent and magnetic descrip¬ 
tion of these complex objects. Asteroseismology has significantly contributed to this 
progress, along with spectropolarimetry and ever improving numerical simulations and 
models of stellar convection and dynamo action. Still much remains to be done. Many 
of the assumptions have been deduced by extrapolating our (incomplete) understand¬ 
ing of our star, the Sun, to other stars. The continuous back and forth comparison 
between the Sun and other stars has been most productive and will continue being so 
with many projects relying both on ground-based telescopes and space-b orn missions 
in the planning.The launch of ESA’s M3 Cosmic Vision PLATO in 2024 ( R.auer et al.l 
l2013l i will continue the rich tradition of data collection initiated by its predecessors 
Corot and Kepler, helping refine our view of stellar dynamics and magnetism. PLATO 
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has a field of view of more than 2200 square degrees and will observe two different 
fields for two years. It will combine huge statistics (about 300,000 stars) with high 
temporal resolution and length of observations to properly extract internal and surface 
stellar rotations as well as signatures of the existence of stellar cycles. Because one 
of the PLATO fields will include the Kepler field, long-term magnetic variability (of 
the order of decades) will be potentially unveiled. Fields of shorter observing times of 
six months (covering in total half of the sky) will add around a million new stars for 
which oscillations will be measured. This large amount of high-quality seismic data 
will be a treasure trove for potential new discoveries and will put new constraints on 
the theories of stellar activity. In the more immediate future, in 2017, the launch of 
NASA’s new planet hunting mission TESS is planned. Although this mission is much 
smaller in scope than Kepler, TESS will observe all stars in the sky down to magnitude 
12.5. Apart from the polar regions of the sky (observed for six months), other fields - 
where the density of stars is higher - will be observed for around 30 days. TESS will 
therefore contribute significantly to stellar rotation period studies, but it is by design 
focused on shorter timescale variability studies than PLATO. 

Furthermore, other recently launched space missions have already or will soon pro¬ 
vide asterosesimic measurements. Apart from the Canadian mission MOST, whose typ¬ 
ical 3 weeks observations per target are not necessarily suited for the type of studies we 
have discussed here, the BRITE constellation of six nano-satellites (http://www.brite-constellation.at) 
and Kepler’s K2 missions are also contributing to the generation new constraints. The 
latter will be particularly important because it provides the ability to employ the new 
tools of asteroseismology in well-studied star clusters, such as M67, the Pleiades, and 
the Hyades. K2 will also sample a variety of stellar populations in the Galaxy, impor¬ 
tant for our understanding of stellar population effects. 

In addition to the space-borne ins truments, the Danish- led ground-based Stellar 
Observations Network Group fSONG. lOrundahl et al.l ll2014l ~)i will soon provide very 
high-quality spectroscopic asterosesismic data of distant stars. Asterosesismic data ob¬ 
served in velocity posses a much higher signal to noise ratio compared to the photo- 
moetric data obtained in space, but for much less stars though. Two of the planned 
eight observing sites, distributed all over the Earth, are already funded, and either in 
construction (China) or already in operation (Teneriffe). The very high-quality data 
expected from SONG will allow us to obtain, for example, information from the very 
outer stellar layers with unprecedented accuracy, information which is much needed to 
improve our knowledge about the still ill-understood surface-effects (see Section 1.2) 
and the details of the physical mechanisms that drive solar-like oscillations. 

The detailed analysis of the legacy star catalog observed by Kepler, the exciting and 
newly planned asteroseismic observation campaigns, the improvements in computation 
power and, last but not the least, our own improved understanding of the physical 
processes in stellar interior will not only bring a wealth of new discoveries, but will also 
bring stellar physics to a new level of sophistication. 

On the theoretical side, recent advances based on analytical work and multi-D 
models of stars of various mass, rotation rates and age, have started paving the way to 
a direct comparison with asteroseismic and spectropolarimetric data. Simulations cou¬ 
pling the angular evolution and magnetic activity on long time scale are in the making 
and should help explain the rotation history of stars. Acoustic and gravity modes ex¬ 
citation and propagation are now computed in detail and stellar cycles have now been 
found in many self-consistent convective dynamo simulations. General trends such as 
getting higher activity levels and shorter cycle period for faster rotating stars have been 
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recovered. Still the exact dynamo mechanisms at work and resulting magnetic prop¬ 
erties (topology, intensity, variability) for any given star are not yet fnlly elncidated. 
This comes about because the relative importance, spatial and temporal phasing and 
nonlinear interactions of the identified processes (e.g. shear in the tachocline, a and 
w effects, stability of the magnetic wreaths, role of the meridional flow and turbulent 
pumping) at the origin of the organisation of the global (large scale) dynamo and the 
establishment of a cyclic magnetic activity vary. Thus the characterisation of the ex¬ 
act dependence of the stellar dynamo properties to stellar parameters such as mass, 
rotation/age and stratification is still an evolving understanding. In order to progress 
further theoretical models and high performance multi-D numerical simulations will 
benefit from having the following inputs/constraints: 

— better constrained variations of the differential rotation in solar-like stars as a 
function of stellar parameters (mass, rotation/age, metallicity) and its latitudinal 
and radial profiles 

— improved characterization of the Pcyc vs Prot relation 

— some assesment of the surface horizontal flows 

— convective power and extent of the convective envelopes 

— link between star spot numbers and stellar parameters, relative contribution to 
stellar brightness between dark spot, and bright faculae 

— high cadence spectropolarimetric maps over the course of many years hence possibly 
mapping stellar cycles for many different stars 

— calibration with asteroseismology of gyrochronology by studying stellar clusters of 
various ages 

It is our hope that in the next decade or so, these multitude of information will 
become available and will aid the community in constraining and further developing 
theories and computational models of solar and stellar activity. 
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